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We review the various methods which have been employed recently to de- 
scribe the thermodynamics of the high temperature quark-gluon plasma 
using weak coupling techniques, and we compare their results with those 
of most recent lattice gauge calculations. Many of the difficulties encoun- 
tered with perturbation theory at finite temperature are in fact not spe- 
cific to QCD but are present in any field theory at finite temperature and 
will be discussed first in the simple example of the scalar field theory. 
We discuss the merits and limitations of various techniques which have 
been used to go beyond perturbation theory in the soft sector, such as di- 
mensional reduction, screened perturbation theory or hard-thermal-loop 
perturbation theory, and $— derivable approximations. In the last part of 
the review, we focus on the later, which lead to a remarkably simple ex- 
pression for the entropy of the quark-gluon plasma. When complemented 
with further, physically motivated, approximations, this approach repro- 
duces accurately the entropy obtained from lattice gauge calculations at 
temperatures above 2.5T C , where T c is the deconfinement temperature. 
This calculation thus provides also support to the physical picture of the 
quark-gluon plasma as a gas of weakly interacting quasiparticles. 
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1. Introduction 

Much of the motivation for studying heavy-ion collisions at ultra-relativistic 
energies is based on the expectation that matter at high temperature and/or 
density becomes simple: because of the asymptotic freedom of Quantum 
Chromodynamics (QCD), one expects that in this regime matter becomes 
a weakly interacting gas of quarks and gluons whose properties should be 
calculable in a weak coupling expansion. 

The existence of weakly interacting quark matter was indeed anticipated 



theoretical evidences for the existence of the quark-gluon plasma are coming 
from lattice gauge calculations (for recent reviews see e.g. RefsPEJ). These 
are at present the unique tools allowing a detailed study of the transition 
region where various interesting phenomena are taking place, such as colour 
deconfinement or chiral symmetry restoration. In this review, we shall not 
consider this transition region, a , but focus rather on the high temperature 
phase, exploiting the fact that at sufficiently high temperature the effective 
gauge coupling constant should be small enough to allow for weak coupling 
calculationjn*™. 

The physical picture which emerges from these weak coupling calcu- 
lations is a simple one, and in many respect the quark-gluon plasma is 
very much like an ordinary electromagnetic plasma in the ultrarelativis- 
tic regim e (with however specific effects related to the nonabelian gauge 
symmetry '™ ' ' I |) . To zeroth order in an expansion in powers of the cou- 
pling g, the quark gluon plasma is a gas of noninteracting quarks and gluons, 
with typical momenta fc ~ T. This is the ideal quark-gluon plasma. The 
interactions appear to alter only slightly this simple physical picture: they 
turn plasma particles into massive quasiparticles, and generate collective 
modes at small momenta. However, while the relevant degrees of freedom 
have been identifi ed long ago, and effective theories to describe their dy- 
namics are knowrP^EIEEIl the calculation of the thermodynamics of 
the quark-gluon plasma using weak coupling techniques has remained un- 
successful until recently. 

Strict expansions in power of the coupling constant have been pushed 
up to order g 5 , both in QCD and in scalar field theories. In both cases, one 
observes a rather poor apparent convergence, the successive contributions 
oscillating wildly, unless the coupling is very small and the deviation from 



on the basis of asymptotic freedom 




most compelling 



For new ideas on how to cover this region by effective field theory methods see Refs. 
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ideal-gas behaviour correspondingly minute. As we shall see, the difficulty 
may be related to the fact that what determines the accuracy of a weak 
coupling calculation in thermal field theories is not only the strength of the 
coupling, but also the magnitude of the thermal fluctuations. These vary 
according to the relevant momentum scales, so that the accuracy of the 
weak coupling expansion depends on which momentum scale contributes 
dominantly to the quantity under consideration. 

To get orientation into these effects, let us then consider a massless 
scalar field theory with g 2 (f> 4 interactions. In the noninteracting case the 
magnitude of the thermal fluctuations of the field is given by the simple 
formula (with Ek = k): 



where n(sk) = l/(e^ £k — 1) is the Bose-Einstein distribution function. We 
shall use this formula as an approximation for the magnitude of the fluctu- 
ations also in the interacting case, eventually modifying appropriately. 
In order to get a criterion for the expected validity of perturbation theory, 
we shall compare the expectation value of the "kinetic energy" (($</>) 2 ) with 
the "interaction energy" which we approximate as g 2 (<f> 2 ) 2 ■ 

For the plasma particles = k ~ T and (4> 2 )t ~ T 2 . Then indeed, when 
g is small, the kinetic energy of a typical particle, k 2 ~ T 2 , is large com- 
pared to its interaction energy ~ y{<f> )t- The short wavelength, or hard, 
fluctuations produce a small perturbation on the motion of a plasma parti- 
cle. However, this is not so for an excitation at the momentum scale k ~ gT: 
then the kinetic energy k 2 ~ {gT) 2 is comparable to the contribution to the 
interaction energy coming from the coupling to the hard modes, g {<fi 
Thus, the properties of an excitation with momentum gT are expected to 
be nonperturbatively renormalized by the hard thermal fluctuations. And 
indeed, the scale gT is that at which collective phenomena develop. The 
emergence of the Debye screening mass ~ gT in a plasma is one of the 
simplest examples of such phenomena. More generally, this renormalization 
of the soft mode dynamics by hard mode contributions invites a description 
in terms of an effective theory for the soft modes, the parameters of this 
effective theory being determined by the hard modes. The building blocks 
of such an effective theory have been dubbed hard thermal loops (HTL). In 
the scalar field theory, there is only one such HTL, which is a mass. But in 
QCD, there is an infinite number of therrM^. 

Let us now consider the fluctuations at the soft scale gT <C T. These 




(1) 
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fluctuations can be accurately described by classical fields. In fact the asso- 
ciated occupation numbers n(sk) are large, and accordingly one can replace 
n(efc) by T/et in eq. Introducing an upper cut-off gT in the momentum 
integral, one then gets: 



Thus g 2 {4> 2 ) g T ~ g 3 T 2 is still of higher order than the kinetic term g 2 T 2 . 
In that sense the soft modes with k ~ gT are still perturbative. Note 
however that they generate contributions to physical observables which are 
not analytic in g 2 , as shown by the example of the order g 3 contribution to 
the energy density, which goes schematically as follows: 



In fact, in contrast to the perturbation theory for the hard modes, the 
parameter which controls the perturbative expansion for the soft modes is 
g instead of g 2 so that is it less accurate in general. As we shall see, much of 
the difficulties with the strict perturbative expansion can be attributed to 
the poor accuracy of perturbation theory in the soft sector, and the present 
review describes various solutions proposed to get around this problem. 

The previous arguments concerning the fluctuations can be extended to 
QCD, where one should compare the kinetic energy (dA) 2 with the magni- 
tude of the nonabelian interactions g 2 A 4 (with A the gauge potential). For 
the plasma particles = k ~ T and (A 2 )t ~ T 2 . The associated electric 
(or magnetic) field fluctuations are (E 2 )t ~ {{dA) 2 )T ~ k 2 (A 2 )r ~ T 4 and 
are a dominant contribution to the plasma energy density. As already men- 
tioned, these short wavelength, or hard, gauge field fluctuations produce 
a small perturbation on the motion of a plasma particle. However, this is 
not so for an excitation at the momentum scale k ~ gT which is non- 
perturbatively renormalized by the hard thermal fluctuations. And again 
perturbation theory for the soft modes is governed by g, not g 2 . 

A new feature emerges however in QCD: at the ultrasoft momentum 
scale, k ~ g 2 T, the unscreened magnetic fluctuations play a dominant 
role. At that scale it becomes necessary to distinguish the electric and 
the magnetic sectors (which provide comparable contributions at the scale 
gT). The electric fluctuations decouple because the Debye screening mass 
m 2 D ^> k 2 ~ {g 2 T) 2 and their contribution is negligible, of order g A T 2 . 
However, because of the absence of static screening in the magnetic sector, 
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we have there ~ k and 



(A 2 ) g2T ~ T 




(4) 



so that (d A) 2 ~ g 6 T 4 ~ <? 2 (A 2 ) 2 : the fluctuations are no longer perturba- 
tive. This is the origin of the well-known breakdown of perturbation theory 



Whereas one may argue that these non-perturbative magnetic fluctu- 
ations can contribute to the deviation of the pressure from the ideal gas 
result, as observed on the lattice, this contribution, being of order <? 6 , is 
presumably numerically rather small sufficiently above the transition tem- 
perature. At any rate, it is clear that the magnetic fluctuations have a priori 
no responsibility for the lack of convergence of the weak coupling expansion 
at lower orders (i.e., up to order <? 5 ). This is further confirmed by the fact 
that a similar lack of convergence is seen also in other theories, like the scalar 
field theory. These magnetic fluctuations require genuine non-perturbative 
treatment, such as tha t ba sed on a combination of dimensional reduction 
and lattice calculation^^!. This will be briefly discussed. But most of the 
review will focus on purely analytical resummation techniques which have 
been proposed to improve perturbation theory in the soft sector. 

Among those, the simplest one which attempts to account nonpertur- 
batively for the dominant effect of interactions in the soft sector, which in 
the case of the scalar field is to generat e a mass for the excitations, is the 
so-called "screened perturbation theory' ^ 8 1 1 xhis consists in a reorga- 
nization of the perturbative expansion, based on the following rewriting of 
the Lagrangian: 



with C' Q = Co — (l/2)m 2 </> 2 . A perturbative expansion in terms of screened 
propagators (that is keeping the screening mass to as a parameter, i.e. not 
as a perturbative correction to be expanded out) has been shown to be quite 
stable with good convergence properties. However, to depends on the tem- 
perature, which makes the propagators explicitly temperature dependent, 
and the ultraviolet divergences which occur in high order calculations be- 
come temperature dependent. Thus, at any finite loop order, the ultraviolet 
renormalization gets artificially modified. 

In the case of gauge theory, the effect of the interactions is more com- 
plicated than just generating a mass. But we know how to determine the 
dominant corrections to the self-energies. When the momenta are soft, these 



at order g 6 . 



C = C - -m 2 4> 2 + -m 2 4> 2 + C mt = C + C' r 



(5) 
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are given by the hard thermal loopd^E3 By adding these corrections to the 
tree level Lagrangian, and subtracting them from the interaction part, one 
generates the so-called hard thermal loop perturbation theory (HTLPT 

C = Cq + LhTL — £-HTL + £-int — C Q + C' int . (6) 

The resulting perturbative expansion is made complicated however by the 
nonlocal nature of the hard thermal loop action, and by the necessity of 
introducing temperature dependent counter terms. Also, in such a scheme, 
one is led to use the hard thermal loop approximation in kincmatical 
regimes where it is no longer a justifiable approximation. 

^-derivable approximations^" provide a more natural approach to prop- 
agator renormalization than screened perturbation theory, although the 
general goal is similar. It allows us to exploit fully the spectral information 
that one has about the plasma excit ations. In particular the hard ther- 
mal loops enter as essential ingredientsj^^^ but in contrast to HTLPT, 
they contribute in an essential way only for external momenta where they 
are accurate. One of the main advantages of such approximations is t hat 
they lead to remarkable simplifications in the calculation of the entropy.ESI 
This allows one, in particular, to bypass some of the difficulties of HTLPT. 
Besides, the entropy makes transparent the underlying physical picture of 
quasiparticles, which (in simpler models) has been shown to provide rather 
accurate fits to lattice 

datJSZEHESEDl. 

The quasiparticle picture assumes 
that the dominant effect of the interactions can be incorporated in the 
spectral properties of suitably defined quasiparticles with small residual 
interactions; this is precisely what comes out of the entropy calculation. 

Finally, we should mention other attempts to improve the behaviour of 
perturbation theory by using various mathematical extrapolation schemes. 
Thus, extrapolations have been constructed, b ased on t he first terms of 
the perturbati ve s eries, using Pade approximant or Borel summa- 
tion techniques^! The resulting expressions are indeed smooth functions 
of the coupling, better behaved than polynomial approximations truncated 
at order g 5 or lower, with a weak dependence on the renormalization scale. 
However these techniques, which offer little physical insight, will not be 
discussed further here. 

The outline of the paper is as follows. In the next section, we shall review 
various features of perturbation theory in quantum field theory at finite 
temperature, using the simple example of the scalar field. This will offer us 
the possibility to illustrate most of the difficulties encountered in thermal 
gauge theories, and also the various resummation schemes which have been 
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proposed in the literature. Section 3 presents a brief summary of the state 
of the art concerning the perturbative calculations in QCD, those based on 
dimensional reduction, as well as recent lattice data. In section 4 we describe 
the calculation of the entropy and emphasize its technical simplicity and 
its physical content. The last section of the review summarizes the various 
conclusions and puts the different approaches into a general perspective. 

2. The scalar field theory as a pedagogical example 

As mentioned in the introduction, many of the difficulties of weak coupling 
calculations at high temperature are not specific to gauge theories. They 
will be illustrated in this section by studying the thermodynamics of a 
scalar field with Lagrangian 

C= 1 -d^d^-^-^. (7) 

We shall often use g 2 = A/4! as an alternative notation for the coupling 
strength. 

2.1. Perturbation theory and its difficulties 

The weak-coupling expansion of the pres sure in th e case m = has been 
computed to order A 5 / 2 (or g 5 ), and readJ^EHEl 



P = Po 



! 15 /g\2 i 15 (g\3 , 135 p , n AnAo \ fg\* 



ttJ 2 



, (8) 



where Pq — (tt 2 /90)T 4 is the pressure of an ideal gas of free massless bosons, 
and <7 2 (/l) = \{fl)/2A is the MS coupling constant at the renormalization 
scale p,. This formula calls for several remarks. 

First we note that, to order A 5 / 2 , P is formally independent of p, (in the 
sense that dP/dp involves terms of order A 3 at least). This is easily veri- 
fied using the renormalization group equation satisfied by the renormalized 
coupling X(p): 

But the approximate expression (JHJ for P does depend numerically on p. 
In order to avoid large logarithms ~ ln(/2/27rT) in the expansion (JHJ one 
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g(p = 27rT) 

Fig. 1. Weak-coupling expansion through orders g 2 , g 3 , g 4 , and g 5 for the pressure 
normalized to that of an ideal g function of g{2itT) in </> 4 theory. 

may choose ft ~ 2ttT. With this choice A (ft) becomes effectively a function 
of the temperature. In the case of QCD where the /3-function is negative, 
this leads to the expectation that the effective coupling becomes small at 
large temperature. 

In Fig. ^ we show the successive perturbative approximations to P/ Pq 
as a function of g(2nT). Each partial sum is shown as a band obtained 
by varying ft from irT to 47rT, which gives an indication of the theoretical 
uncertainty in the calculation. Thus, in drawing this plot, g(fi) has been 
calculated from g(2irT) by using the renormalization group equation J5|J|. 
The lack of convergence of the weak-coupling expansion is evident in Fig.^ 
We can infer from Fig. ^that the expansion in powers of g ceases to make 
sense as soon as g(2nT) > 0.4. 

Another remark concerns the presence of terms non-analytic in the 
coupling constant A, i.e., terms of order g 3 ~ A 3 / 2 , g 5 ~ A 5 / 2 , and also 
terms involving In g, in the perturbative expansion in eq. (JSJ. As we shall 
see shortly, these terms result from infinite resummations of subsets of 
diagrams of perturbation theory. 

At this stage, it is instructive to compare with the corresponding ex- 
pansion of the screening "Debye" mass mp, defined from the the static 
propagator by b 

p 2 + m 2 D + 11(0, p) = at p 2 = -m 2 D , (10) 



b A definition through the location of the pole rather than the infrared lim it o f the 
self-energy is crucial to make mrj gauge- indep endent in non-Abelian theorieJSHj this 
definition is also renormalization-group invarianJ22E21 
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m/T 1 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 

g(p = 27rT) 

Fig. 2. Weak-coupling expansion to orders g 2 , g 3 , and g 4 for the screening mass nor- 
malized to the temperature as a function of g(2irT). 

where H(po, p) is the self-energy. For the massless case (m = 0), the weak 
coupling expansion of m 2 D is known to order g 4 , and r eadJ23 

m 2 D =g 2 T 2 {l-Z 9 -^- 

V IT l 

Note that, with our notation g 2 = A/24 for the coupling constant, the 
leading-order contribution in eq. is simply 

m D =gT, (12) 

which happens to be the same relation as in pure-glue QCD to leading 
order, cf. I|91l) . As for the pressure, non-analytic terms (g 3 , In 1/g) occur in 
the expression of m 2 D . 

In Fig. El we show the screening mass normalized to the leading 
order result rhjj as a function of g(2irT), for each of the three successive 
approximations to m 2 D . As in Fig. Q the bands correspond to varying /2 
from 7rT to AnT . The poor convergence is again evident, with a pattern 
similar to that in Fig.Q] a large difference between the order-c/ 2 and order-g 3 
approximations, and a larger sensitivity to the value of the renormalization 
scale for the order g 4 . 

In order to get some insight on this peculiar behaviour of the expansions 
© and HI If) . we shall briefly recall how the various terms are obtained, up 
to order g 3 . 

We start with the zeroth order calculation. We have Pq = — CIq/V, with 

1 rp f d 3 k , / 2 ,i 2 j 2n 
-y^2 T 2^ 7CT ln K + k +m) 

n v J 




In — In - - 2.415 

27rT 3 7T 



O 



(3 3 )}-(ll) 
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-\J&« + T J&"*--~ IT >- (13 > 

where oj n = 2nirT is a bosonic Matsubara frequency and e\ = k 2 + m 2 . The 
first term in the second line of eq. I|13|l is the sum of the 'zero-point' energies 
of free massive bosons. This is ultraviolet divergent, but independent of 
the temperature; it can therefore be absorbed in the redefinition of the 
vacuum energy, and discarded. The second term is temperature-dependent 
and finite; this is recognized as minus the pressure of the ideal gas of scalar 
particles. For m <§; T this can be expanded as: 

f d 3 k , , Bf . tt 2 T 4 m 2 T 2 m 3 T . . 

where the neglected terms start at order m i ln(m/T). 

Consider now the contributions of order g 2 to the self-energy and the 
thermodynamic potential, denoted as 1T2 and O2, respectively. 1T2 is given 
by the "tadpole" diagram of Fig. [3] This gives: 

Ho = — T / - — — = 7: x = — I(m) 

2 ^ J (2tt) 3 uj 2 n + k 2 + m 2 2 V ' 

A f d 3 k l + 2n(e k ) _ A 
= 2 7 (2^ 2 £fc = 2 [/o(m) + /T(m)] ' (15) 

where n(o;) is the Bose-Einstein thermal distribution: 

n M = c p^_ 1 > ( 16 ) 

and the notation 7(m) for the sum- integral, as well as Io(m) and It(iti) 
for its zero temperature and thermal contributions, will be used repeatedly 
in the following. The first term in the second line of eq. (|15|1 is ultraviolet 
divergent, but independent of the temperature, so it can be absorbed into 
the definition of the (zero-temperature) mass. For instance, is we choose 
the mass parameter m in the Lagrangian to be the physical mass at T = 0, 
we must add a counterterm ^8m 2 (j) 2 in the Lagrangian, with Sm 2 chosen 




Fig. 3. Lowest order correction to the self-energy in scalar 4 theory. 
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^2 = ^lT( m )^-\-- — -—- 2 [ln 7J -^ + 2 1 -l)l (18) 



Fig. 4. The lowest order correction to the thermodynamic potential in cj> 4 theory. 

so that: 

^I (m) + 8m 2 = 0. (17) 

Then, the lowest-order self-energy correction due to the thermal fluctua- 
tions reads: 

A . A f T 2 mT 
2 y ' 2 \ 12 4vr 16tt 2 V (4ttT) 2 

where the expansion in the r.h.s. holds when m < T. In particular, when 
m = 0, we recover the result IT2 = m 2 D , cf. eqs. (|11|I - H12|) . This quantity is 
generally dubbed a "hard thermal loop" (HTL), because when m the 
dominant momenta in the one-loop integral (|18fl are hard (k ~ T). 

The contribution O2 is given by the 2-loop diagram in Fig. 01 and reads: 

T = + \ 6m2{4,2)o = 5 (/|(w) " 7 ° 2(m)) ' (19) 

where we have also included the contribution of the mass counterterm and 
used eq. (|17|) to verify that the temperature-dependent ultraviolet diver- 
gences cancel out, as they should. The second term in the last line in eq. (|19(1 
is divergent, but contributes only to the vacuum energy density, so that it 
can be discarded. For m = 0, eq. (|19|) yields: 

*~iis r '-T(?)"*- « 

in agreement with eq. (Hj. 

At the next order in the loop expansion, the self-energy receives con- 
tributions from the 2-loop diagrams. However, some of these contributions 
are infrared divergent in the massless limit m — 0. Consider, e.g., the second 
diagram displayed in Fig.0 The upper loop (together with the correspond- 
ing mass counterterm, not shown explicitly in Fig. EJ) is recognized as the 
"hard thermal loop" computed previously. Thus, this particular 2-loop con- 
tribution reads: 

nW - A rV f d3fc (21) 
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The three-dimensional integral over k in eq. I|21|) has a linear infrared 
divergence, coming from the term with uj n = 0. (For the other terms, 
uj n 7^ 0, the Matsubara frequency w n acts effectively as an infrared cutoff.) 
A non-zero mass to would cut off this divergence, and give a contribution 
~ Xihjj (T/m) ~ g 4 T 3 /m. As we shall see, when to <C gT, the real cut- 
off is not actually the mass m itself, but the "thermal mass" rhry = gT 
induced by the thermal fluctuations. Then IT^ ~ g A (T 3 /m D ) ~ g 3 T 2 . 
This illustrates how nonanalytic terms arise in eqs. JHJ and potential 
infrared divergences are cut-off by 'thermal masses' which, through their 
dependence on g, modify the original expansion in powers of the coupling 
constant. We discuss this in more detail in the next subsection. 

2.2. Order g 3 from various resummation schemes 
2.2.1. The next-to-leading-order thermal mass 

The effect of a thermal mass is most directly seen by carrying out pertur- 
bative calculations with the dressed propagator 

For the calculation of the tadpole, this corresponds to resumming the di- 
agrams of bare perturbation theory shown in Fig. [5] (each of the inserted 
tadpoles in Fig. [5] is accompanied by a mass counterterm which removes its 
ultraviolet divergence in the vacuum, cf. eq. (|17|l ). 

Fig. 5. Diagrams that are resummed when computing the tadpole in Fig. l3l with the 
dressed propagator in eq. 1221 . 

Since m/j = gT, this resummation involves arbitrarily high powers of 
g. However, the final result can be read from eq. I|18|l in which we set 
to = to ry = gT: 

ITadpoic = g 2 T 2 - -g 3 T 2 + 0{g 4 \ng) . (23) 

This formula exhibits the next-to-leading order contribution to the thermal 
mass, of order g 3 , consistent with eq. lflT|l . This order-g 3 contribution comes 
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from the term linear in m in the r.h.s. of eq. (JTSJ. It thus is related to the 
non-analyticity of the thermal integral i-r(m) as a function of m 2 and to 
the infrared divergences which occur when one attempts to evaluate this 
integral after expanding the integrand in powers of m 2 . These infrared 
divergences are easily identified in the sum of the diagrams in Fig. El with 
at least two loops: 

A ^ f d 3 k f 1 1 

TT ■ — — T x 



2 (2tt) 3 \ul + k 2 + m 2 D ul + k 2 

= 12 9 2 TV 1j W-=%V , (24) 

where the subtracted term is the massless tadpole. The term with I = 1 
corresponds to the 2-loop diagram that we have considered separately in 
eq. (J2J. Each of the terms with uj 2 > in eq. (|24(l starts to contribute 
at order g 2 m 2 D ~ g 4 . But the terms with uj n = are individually infrared 
divergent for any /, and the divergences cancel only in the sum of all the 
terms. In fact, the contribution of order g 3 arises from the static modes 
alone, i.e., the modes with zero Matsubara frequency. If we separate their 
contribution in eq. I|24|) . we obtain indeed 

^- J wivhs-i?} (25 » 

where we recognize the order g 3 contribution in eq. I|23|l . 

It is instructive to look at this contribution from another perspective. 
As we have argued before, the dominant contribution to 7r(?n) in eq. (|18(l 
comes from momenta k ~ T. The subleading contribution is coming from 
soft momenta k -C T. To isolate it, we subtract the hard modes contribu- 
tion, which amounts to subtract ir(0) from /t(ti) in eq. . and write: 

After multiplication by A/2 = 12<? 2 , one recognizes the g 3 term of eq. I|23(l . 
In writing the approximate equality above, we have used the fact that the 
dominant contribution in the integral (|26|l comes from momenta k -C T, so 
that we can use 

n{e k ) S3 — for £ fc <T. (27) 

£k 

Note that the above resummation of the thermal mass, focusing on 
the static mode, hides an ultraviolet problem that we have ignored so far. 
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Namely, in order to obtain the finite result in eq. (|23[) . we had to cancel an 
ultraviolet divergence, cf. eq. (|17|l . which now involves the thermal mass, 
and is therefore temperature dependent. The T-dependent divergent piece 
is of order g 4 however (since proportional to g 2 m|,), and thus does not 
formally affect the order-g 3 calculation. But it signals that our calculation 
is not complete. The divergence is identified in eq.J2H) in the term with 
I = 1, i.e., in the 2-loop diagram in Fig.0 after performing the Matsubara 
sum. The overall divergence of the 2-loop diagram arises as the product 
of a vacuum divergence in the 4-point function subgraph (the lower loop) 
times the finite-temperature piece of the upper loop. Thus, the divergence 
is removed by the vacuum renormalization of the 4-point function, i.e., by 
adding a 1-loop counterterm to the vertex in the tadpole. 

2.2.2. The plasmon effect in the pressure 





Fig. 6. The "ring" diagrams responsible for the order-g 3 effect in the pressure. 

The contribution of order g 3 to the pressure, also known as the "plasmon 
effect" , comes again from resumming diagrams which are infrared divergent 
in the naive diagrammatic expansion. These are the "ring" diagrams in 
Fig. in which the central loop is static (u> n = 0), while each outer loop 
represents an insertion of rh 2 D . This series can be easily summed up, with 
an infrared finite result 



T 
2" 

T 
1> 



d 3 k 
(2^)3 

d 3 k 
(2tt) 3 



E (-D 

l>2 



l-l 



' D 



k 2 



rh 3 D T 
12tt 



which is the order-g 3 contribution in eq. 



(28) 



2.2.3. Static resummation 

From the previous examples, we learned that (a) the occurrence of infrared 
divergences in the diagrammatic expansion calls for resummations and (b) 
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it appears sufficient, in order to cure these divergences, to resum rh 2 D in the 
static Matsubara modes alone. 

Since m 2 D summarizes interaction effects, care should be taken to avoid 
overcounting in carrying out high order calculations within that scheme. 
One way to proceed systematically is to add and subtract a mass term for 
the static modes in the (imaginary-time) Lagrangian: 

£^i(^) 2 + i(V0) 2 + A 4 

= \{dA? + \{Vcp) 2 + \ + ±<f>* + A£j (29) 

where AC = — ^w 2 ^ 2 , an d ^o( x ) = T J Q dr0(r, x). This allows for a re- 
organization of the perturbation theory in which the tree-level amplitudes 
are generated by the first three terms in eq. I|29[l while AC is treated for- 
mally as an additional interaction, of order g 2 , to ensure that the HTL is 
not double counted in the static sector. 

This procedure preserves a strict correspondence between diagrams 
(or well identified pieces of it) and powers of g. Thus, it is well suited for 
weak coupling expansions of static quantities, and it is indeed the scheme 
in which the higher-order calculations for the pressure and the screen- 
ing mass alluded to befo re have b een originally performed (for both scalar 
theoryEMSIIand QCE0D22E3I). 

2.2.4. Dimensional reduction 

The special role of the static mode 0o may be handled in a more systematic 
fashion by constructing an effective theory in which the non-static modes 
4>v=£o are 'integrated out' in perturbation theory. Specifically, one can write 
the thermal partition function as the following path-integral: 

Z=N y"p(0o)exp{-5 cff [0 o ]} , (30) 

where 

cxp{-S cff [0 o ]} = M' J 2?(</v /0 )exp j- J\t J tfxL B {x)\ , (31) 



and 5 ff [4>o] is the effective action. Aside from the direct classical field con- 
tribution (i.e., the restriction of eq. JJJ) to the static mode (f>o), this effective 



c This is true, strictly speaking, only if one uses dimensional regularization to control the 
ultraviolet divergences. 
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action receives also contributions which, diagrammatically, correspond to 
connected diagrams whose external lines are associated to <fro, and the inter- 
nal lines are the propagators of the non-static modes 4> u ^o- Thus, a priori, 
S c s[<Po] contains operators of arbitrarily high order in <fr , which are also 
non-local. In practice, however, one wishes to expand S B s[4>o] in terms of 
local operators, i.e., operators with the schematic structure a m) n V m 0Q with 
coefficients a m , n to be computed in perturbation theory. 

It is then useful to introduce an intermediate scale A (A < T) which 
separates hard (k > A) and soft (k < A) momenta. All the non-static modes, 
as well as the static ones with k > A are hard (since K 2 = ui^ + k 2 > A 2 
for these modes), while the static (oj n — 0) modes with k < A are soft. 
Thus, strictly speaking, in the construction of the effective theory along the 
lines indicated above, one has to integrate out also the static modes with 
k > A. The benefits of this separation of scales are that (a) the resulting 
effective action for the soft fields can be made local (since the initially non- 
local amplitudes can be expanded out in powers of p/K, where p <C A 
is a typical external momentum, and K > A is a hard momentum on 
an internal line), and (6) the effective theory is now used exclusively at 
soft momenta. This strategy, which consists of integrating out the non- 
static modes in perturbation theory in order to obtain an effective three- 
dimensional theory for the soft static mod es, is generally referred to as 
"dimensional rediictio n '™ ™BB™l 

As an illustration let us consider a massless scalar theory with quartic 
interactions; that is, eq. JjJ with m — 0. The effective action for the soft 
fields takes the form: 

SM = ^(A) + J d 3 x ji(V<M 2 + \m 2 (A)^ 2 + g 2 (A)4 + •••}, (32) 

where J~{A) is the contribution of the hard modes to the free-energy, and 
the dots stand for higher-order, local operators. We have changed the nor- 
malization of the field (</>o — > VTcpo) so as to absorb the factor /3 in front of 
the effective action. The parameters of the effective Lagrangian in eq. I|32[l , 
like M 2 (A), 33(A), etc., are computed in perturbation theory (they involve 
diagrams in which all the internal lines are hard), and depend on the sep- 
aration scale A, the original coupling g, and the temperature T (note that 
the effective theory depends on the temperature only through these pa- 
rameters). To lowest order in g, g 2 w g 2 T, and M ~ gT, as we shall see 
shortly. 

Note that the scale A acts as an explicit ultraviolet (UV) cutoff for the 
loop integrals in the effective theory. Since it is arbitrary, the dependence on 
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A of soft loop contributions must cancel against the dependence on A of the 
parameters in the effective action. Let us verify this cancellation explicitly 
in the case of the thermal mass M of the scalar field, and to lowest order 
in perturbation theory. To this order, the mass parameter M 2 (A) in the 
effective action is obtained by integrating over hard momenta within the 
one-loop diagram in fig.|3|(cf. eq. j!5|)). This gives 



The first two terms within the last integral are the same as in eq. (|f 5(1 with 
m = 0. The first of them, involving the thermal distribution, gives the HTL 
contribution rh 2 D — g 2 T 2 , while the second one, involving l/2fc, is removed 
by the vacuum renormalization. The third term, involving the ^-function, 
is easily evaluated. One finally gets: 



The A-dependent term above is subleading, by a factor A/T <C 1. 

The one-loop correction to the thermal mass within the effective theory 
is given by the same tadpole diagram, fig. |3| evaluated with the massive 
propagator l/(fc 2 + M 2 (A)), coupling constant g 2 « g 2 T, and ultraviolet 
cut-off A. We obtain 



where the terms neglected in the last step are of higher order in rho /A or 
A/T. 

As anticipated, the A-dependent terms cancel in the sum M 2 = M 2 (A) + 
SM 2 (A), which then provides the physical thermal mass within the present 
accuracy: 




(33) 




(34) 




M 2 = Af 2 (A) + SM 2 



(A) = g 2 T 2 




(36) 



which agrees again with eq. 123J1 . 
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2.2.5. Screened perturbation theory 

Still another strategy has been used for taking into account the effects 
of the screening mass to all orders, rather than perturbatively. It consists 
of resumming the HTL rh 2 D in all the Matsubara modes, static and non- 
static (as we did in subsection 12.2. 1|) . It can be formulated via a simple 
generalization of eq. Q29J1. viz. 

C = -d^t - ~ rh 2 D 4> 2 - + A£; AC = 1 ^2 ) (3?) 

with AC to be treated as a quantity of order g 2 in perturbative calcula- 
tions. This scheme is commonly referred to as screened perturbation theory 
To appreciate its virtues, and its limitations, it is again best 
to work out some examples. 

We have obtained earlier the expression of the tadpole calculated with 
the dressed propagator (f2"2")l (cf. eqs. (|2"3"|) and fH^l 'l: 

Iltadpoic - ^9 2 lT{rh D ). (38) 

We have seen that n ta dpoie(wD = gT) 1 when expanded in powers of <?, re- 
produce the correct perturbative result to order g 3 (see eq. But, as 
g grows, it deviates rapidly from this perturbative result, as can be seen 
in Fig. \7\ The poor convergence of the weak coupling expansion exhib- 
ited here is related to that of the high-temperature expansion in eq. (|18l) 
which, as we have seen, involves in fact non analytical terms. By includ- 
ing the screening mass in the tree-level Lagrangian, and not treating it a 
perturbative quantity of order g, SPT provides a smooth extrapolation to 
large values of g, which represents a definite improvement over perturbation 
theory. 

However, in contrast to the static resummation scheme based on 
eq. H29f> . the calculations within SPT generally involve temperature- 
dependent ultraviolet divergences. Of course, these are compensated by the 
"counterterm" AC, but only in all-order calculations. At any finite-order, 
the compensation holds only to the perturbative accuracy d of the calcula- 
tion. For instance, in order to obtain iJBEJ); one has to introduce the mass 
counterterm Sm 2 = — (X/2)Io(rhi)), which depends on T, via tod. Specifi- 
cally (with dimensional regularization, for definiteness, and p, 2 — Aire r /i 2 ): 

' 3 ~ £l - 1 _ 3g 2 m 2 D (2 ^ fl 2 \ 



d By the "perturbative accuracy" of an all-order result we understand the highest order 
in g that is correctly reproduced by the weak coupling expansion of that result. 
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i/T 




Fig. 7. The one-loop thermal mass in SPT, mo = Ilj^ . e , eq. 1381 . normalized 
to the temperature as a function of g (full line). For comparison, we also present the 
perturbative estimates to order g 1 (long-dashed line) and g 3 (short-dashed line). 



where £k{rho) 2 = k 2 



D ' 



This formula makes explicit the fact that 



the T-dependent divergence is of order 0(g 4 ), and thus lies beyond the 
perturbative accuracy of this one-loop calculation, which is 0(g 3 ). 

As further illustration of the difficulties with SPT, let us compute the 
pressure to 2-loop order, and compare the result with its perturbative ex- 
pansion to order g 3 . 

To one-loop order in the effective theory, the thermodynamic potential 
is given by eq. I|13|) with m — ► rhu. That is: 



no 
v 



2 ^ 

n 

d 3 k 



d 3 k 
(2?r)3 



ln(w 2 



(2tt) 3 

with Po(m) given by eq. I] 14(1 : 



£k{m D ) - P Q (rh D ) 



Po{m D ) 



n 2 T 4 
90 



8 Vtt 



15 

T 



0(g 4 lng) 



(40) 



(41) 



This overincludes the order-g 2 effect by a factor of 2 (compare to eq. ©), 
but includes correctly the term of order g 3 (the fact that the plasmon effect 
is included correctly is somewhat accidental: If we were to resum not just 
the leading-order thermal mass rh 2 D: but also the NLO correction to it, 
Amjj = — (3/ir)g 3 T 2 , one would modify the order-g 3 content of f^o)- The 
first piece in eq. I)4()[l. i.e., the sum of the "zero-point energies", is also UV 
divergent, and the divergent terms depends on T since they depend onmn. 
To eliminate these divergences, we may introduce a temperature-dependent 
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counterterm: 

dn 1 f d 3 k , . 1 f d A k , , l2 l2 A ,, 

— = -2JWf £k{mo) = ^\wf Hk ° + k + mV) ' (42) 

whose dependence upon T starts already at order g 2 . (This dependence can 
be pushed to order g 4 , though, by using dimensional regularization.) This 
lowest order calculation reveals another unsatisfactory feature of SPT: the 
tree-level pressure calculated with the Lagrangian i|37|) fails to reproduce 
the lowest-order effect of the interactions. 

Of course these drawbacks are partially corrected when we move to the 
next order, that is to two-loop order, which we do now. At this level, there 
are three contributions: the two-loop diagram in Fig. 0] and two one-loop 
diagrams involving mass counterterms, namely AC in the Lagrangian 13711 . 
whose role is to correct for overcounting, and the ultraviolet counterterm 
Sin 2 , that has been introduced before, in eq. I|39l) . Altogether: 



9l = - 

V ~ 8 



[I(m D )] 2 + \ (Sm 2 - ml) I{m D ) . (43) 



It is convenient to combine the last counterterm with the tree-level expres- 
sion f2 (written as in the first line of eq. j%Q]| ): 

— — — — limn ) 
V 2 y ' 



[ ^ W + k 2 + m 2 D ) g> . a } . (44) 

Clearly, the main effect of AC is to cancel out the whole contribution of 
order g 2 that was originally present in the tree-level pressure. This includes 
the "anomalous" order-g 2 effect within Po(mr)), and the divergent contri- 
bution of order g 2 that was eliminated by the counterterm 142|) at one-loop 
order. The terms in eq. 143|l can be simplified as follows: 

^ [IK)] 2 + ^ I(m D ) = ± [I T (m D )} 2 - ^ [I (m D )} 2 , (45) 

(cf. eq. H19|) ). The last term is divergent and depends on the temperature, 
via rhr), but this dependence counts only to order g 4 , or higher. The order- 
g 2 is obtained by replacing It {mo) by It(0) in the first term in eq. J3SJ. 
The order g 3 receives contributions both from eq. (|44|l (as the lowest order 
correction beyond the ideal gas result), and from eq. H45[) (as the next- 
to-leading order contribution). In both cases, the integration involves soft 
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momenta only. Thus, we can restrict ourselves to the static mode u> n = 
in eq. (@IJ), and get: 

irZ-gMlafl + gg) - \ (46) 

while in eq. (|45|> we may proceed as follows: 

^ (lT(m D ) f - ± (I T (0) f ~ ^7 T (0) (J T (m D ) - 7 T (0)) 

„ 2 1 f d 3 k / 1 1 \ 

where we have also used eq. (|26[1 . Clearly, after adding the two contributions 
in eqs. (|46|) and l|47|) . we are left with the standard expression for P3, 
eq. (J2SI- 

Now, all the remaining ultraviolet divergences in eqs. I|44|l and (|45|) 
are to be cancelled by appropriate counterterms. One could be satisfied 
by the fact that the T-dependence of these counterterms starts at or- 
der g 4 , which is beyond the perturbative accuracy of the present 2-loop 
calculation (namely, 0(g 3 )). But leaving temperature dependence in UV 
counterterms is somewhat unphysical and would cause problems in non- 
perturbative calculations. Moreover, ambiguities arise from the choice of 
the subtraction scheme: One could e.g. either choose to subtract all the 
T-independent terms of the thermal pressure before identifying rho with 
a thermal mass, or one may only subtract pole terms minimally on the 
grounds that also the terms not explicitly dependent on T become so after 
making rho temperature-dependent. 

Some of these drawbacks may actually be absorbed by the flexibility 
in the choice of the mass parameter fhu in 1|37|1 • Obviously, the results of 
SPT would be independent of this mass if they were obtained in an all 
order calculation. But for a finite-order calculation, the choice matters. 
One may try to optimize perturbation theory by a princi ple of minimal 
sensitivity resulting in a self-consistent gap equation ^ 1 '' I This has been 
shown to reduce greatly (though not completely) the dependence on the 
subtraction procedure'^ (see also below, Fig. EJ. Indeed, in Ref. the 
calculation of the pressure of scalar 4> 4 theory has been extended to 2-loop 
order for the screening mass and 3-loop order for the pressure (meaning a 
perturbative accuracy of order g 5 in both cases), with results showing an 
excellent apparent convergence: i.e., the difference between the 2-loop and 
3-loop results for the pressure remains rather small for all g < 1. 
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2.3. Self -consistent resummation 

2.3.1. Skeleton expansion for thermodynamical potential 

A systematic way to take into account in the calculation of the thermody- 
namics, screening effects, or more gener ally propagator renormalisations, 
is to use the representation ^ * ^ I vv l v 4 1 of the thermodynamic potential 
n = —PV in terms of the full propagator D: 



3Q[D] =-logZ=- Tr log D^ 1 - i Tr ILD + &[D], (48) 

where Tr denotes the trace in configuration space, 8 = 1/T, and $[D] is 
the sum of the two-particle irreducible (2PI) "skeleton" diagrams with no 
external legs: 



$[D] = 1/12^^ +1/8 (^)(^) + l/48 




(49) 



(where in the tp 4 example the first diagram is absent of course). The full 
propagator D is expressed in terms of the proper self-energy II by Dyson's 
equation (Do is the free propagator): 

D = i , (50) 

Do +n 

and II itself, which is the sum of the one-particle irreducible (1PI) diagrams 
with two external lines, is obtained from $[D] by: 

5$[D]/5D = hi. (51) 

By using this relation, one can check that f2[D] is stationary under varia- 
tions of D (at fixed Do) around the physical propagator: 

Sn[D}/5D = for D = (D^ 1 + II) -1 . (52) 

We shall usually refer to eq. (|50() in which II is given in terms of D by 
eq. ifSTJl as a 'gap equation'. 

An explicit expression for f2 is obtained by performing the summations 
over the Matsubara frequencies in eq. I|48|l , using standard contour integra- 
tion techniques. One obtains: 

n/V= f-£^n(uj) (lmlog(-w 2 + fc 2 + n) - ImHD) +T<f>[D]/V (53) 

J (27T) 
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where n(ui) = l/(e^" — l), and the imaginary parts are defined with retarded 
prescription. For instance: 

lmD(oj, k) = lmD{u + ie, k) = ^ ® , (54) 

where p{ui, k) is the spectral function, which enters the following represen- 
tation of the analytic propagator D(u>,k), valid for frequencies lo off the 
real axis: 

f°° dk p{k ,k) 

D(u),k)= - — - . 55) 

J-oo 27r k a - u 

The stationarity condition l|52l) can be used as a variational principle 
to deduce non-perturbative approximations for the physical propagator. A 
convenient approximation scheme of this type is the self-consistent (or 
derivable" ) approximation^! obtained by selecting a class of skeletons in 
<£>[£>] and calculating II from eq. I|51|l. We shall see also that the stationarity 
property l|52[l brings simplifications in the calculation of the entropy. 



2.3.2. A simple model 

As a simple illustration of the general formalism of the previous subsection, 
we consider here the 2-loop approximation to <£>, given by the first skeleton 
in the r.h.s. of cq. (|49l) . c The corresponding skeleton for II is the dressed 
tadpole (see Fig. |3J), and is local. Thus, the self-consistent self-energy is 
simply a mass, and we shall write II = m 2 . (The vacuum mass is set to 
zero.) The gap equation itSlJl reads then: 

m 2 = -I(m). (56) 

In contrast to the previous renormalization of eq. I|38|l (recall eq. H39|) ~). 
it is possible to renormalize eq. I|5t)|) without introducing a thermal counter- 
term. This is done through coupling constant renormalization. To proceed, 
we first rewrite eq. (|56|) as (with the same notations as in eq. (|39|I % 1: 

e , m 2 /2 /I 2 \ f d 3 k n(e k ) , s „ s 

= - Ao ^t +l0 4 +1 ) + A »i(^ir + 0(e) ' (57) 



e The s elf-c onsistent solution to 3-loop order has been worked out using mass expansions 
in Ref . with satisfactory numerical results. However, the negative conclusion therein 
concerning renormalizability of ^-derivable approximations beyond two-loop order has 
since been refuted by Refs. 
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and relate Ao to the renormalized coupling A by: 

1 _ n e 1 
A ~ A^ + 16ttV 



(58) 



Then, the equation takes the following, manifestly finite, form (for e —> 0): 



2 A f d 3 k n(ek) Am 2 



-'^/wif^rp-V' (59) 

One should stress a subtle point in the procedure: expressing the bare cou- 
pling Ao in terms of the renormalized one A according to eq. (|58() is not 
sufficient to render finite the r.h.s. of eq. i|57|l for arbitrary values of m 2 , 
but only for these values which satisfies the finite equation l|59[) . 

Note that the 1-loop correction to the 4-point function that is consid- 
ered here is iterated only in one channel among the three possible ones. 
Accordingly, the /3-function for the renormalized coupling constant l)58[l: 

dA *' (60) 



dlog/i 1671" 2 

(this follows by noticing that Ao is independent of /2 in eq. JSSJl) is only 
one third of the lowest-order perturbative /3-function for this scalar field 
theory (see eq. ©)• This is no actual fault since the running of the coupling 
affects the thermodynamic potential only at order A 2 which is beyond the 
perturbative accuracy of the 2-loop <&-derivable approximation. (In order 
to see the correct one-loop /3-function, the approximation for $ would have 
to be pushed to 3-loop order.) By using eq. iffiO")) . one can check that the 
solution m 2 of eq. (|59|l is independent of fi. 

By using the previous formulae, one can easily compute the pressure^! 

fe m 2 /-_^fc_n(e*) _m^_ 
J (2tt)3 10gU )+ 2 J (2tt)3 2e k + 128^ 2 ' 

which differs from the pressure Po(m) of an ideal gas of massive particles, 
eq. H14J1 . by the last two terms in the r.h.s. 

By construction, the 2-loop <!>-derivable approximation is perturba- 
tively correct to order g 3 , and, indeed, it can be easily verified that the 
weak coupling expansions of eqs. 1)59(1 and l|fil(l include the correct pertur- 
bative effects of order g 2 and g 3 (together with incomplete contributions 
of higher orders). But the complete, self-consistent results, as obtained by 
numerical evaluation in eqs. H59[) and ()61[) . show a much better behaviour 
than the respective perturbative expansions to order g 3 , in the sense of be- 
ing monotonic with g (unlike the perturbative approximants of order g 3 ), 
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Fig. 8. Numerical solution of the mass gap equation 1591 with p, = 2ttT in comparison 
with the perturbative result to order g 2 and g 3 for m (long and medium dashed lines, 
resp.). The short-dashed line just below the full line corresponds to the quadratic gap 
equation 1651 following from the requirement of only soft self-consistency. 

and showing saturation at very large g > 2 (in contrast to the order-g 2 esti- 
mates). This behaviour is illustrated in Figs. 0D and [5] for the self-consistent 
mass (|59fl and pressure (|61|) . respectively. 

Also shown in Figgis the result of SPT at two-loop order when the mass 
parameter is determined by a principle of minimal sensitivity. It turns out 
that for a minimal subtraction of the additional UV divergences of SPT and 
identification of the corresponding regularization mass scale with that used 
in the (MS) renormalization of the coupling constant, the result of SPT 
coincides with that of the 2-loop ^-derivable approximation. However, if 
the additional UV divergences of SPT are subtracted differently, this coin- 
cidence no longer occurs. As an example, in Fig.Elthe grey full line displays 
the SPT result with a "maximal" subtraction of the additional UV diver- 
gences where P(T) is made finite by subtracting P(0) when treating rho as 
temperature-independent (since it is just a parameter in the Lagrangian), 
and only afterwards turning it temperature-dependent by fixing m,u by a 
variational principle. Although (for N — 1 scalar theory) one cannot decide 
how the exact result would look like, the result appears to be less satisfac- 
tory in that it exceeds the free pressure for g > 1.5. The prescription that 
has been used in the application of SPT to QCD is in f act that of a minimal 
subtraction of the additional UV divergences 

We have previously argued that the expected accuracy of perturbation 
theory is a priori different in the soft and the hard sectors. Within the 
present solvable model, it is possible to verify this argument explicitly. To 
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Fig. 9. The 2-loop <t>-derivable pressure 1611 in comparison with the perturbative results 
to order g 2 and g 3 (long and medium dashed lines, resp.). The short-dashed line just 
above the full line corresponds to the pressure evaluated perturbatively for hard modes 
and self-consistently only for the soft modes, eqs. 1621 and 1661 with 1651 . SPT at 2-loop 
order turns out to coincide with the 2-loop ^-derivable given by the full line, when the 
former is subtracted minimally, but gives different results (grey full line) when subtracted 
"maximally" (see text). 

this aim, we shall reformulate the self-consistent calculation of the thermal 
mass and the pressure in an approximate way, exploiting the distinction 
between the hard and soft sectors. That is, we shall build on the strategy 
of effective theories to obtain perturbatively, through a calculation in the 
hard sector, the coefficients of an effective Lagrangian for the soft sector. 
We shall then look for an "exact" solution in the soft sector, which amounts 
here to use the two-loop self-consistent approximation. As we shall see, this 
strategy, which mimics that used in calculations based on dimensional re- 
duction, turns out to be remarkably accurate. It also provides a clearer 
physical justification for approximations that have been motivated previ- 
ously by other considerations. 

Thus, we shall approximate the fully self-consistent 2-loop result for 
the pressure, eq. I|61|) . by P = P( 2 ) + Ps, where P( 2 ) is the perturbative 
result to order g 2 , which is due to the hard modes alone: 



and P s is obtained via a 2-loop self-consistent calculation within the ef- 
fective three-dimensional theory with Lagrangian To the order of 
interest, one can set M 2 ss m 2 D and g% « g 2 T in eq. and neglect all 
the other interaction vertices. 




(62) 
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The corresponding free-energy functional reads then: 



V ZJkl k 2 + m|, + II 

+ K/ fefc2 + 4 + n s ) ■ (63) 

where J, denotes the momentum integral in 3 — d dimensions. Note that, as 
compared to Sect. l2.2~4l here we find it more convenient to use dimensional 
regularization, with the scale fi chosen as the separation scale A between 
soft and hard momenta. This has the advantage that all the power-like 
divergences in integrals like those in eq. I|63[l are now set to zero, so, con- 
versely, the parameters of the effective Lagrangian are independent of A to 
the accuracy of interest 2 -'- 1 . 

In eq. H63[l . n s is the scalar self-energy within the 3 — d effective theory, 
and is determined from a gap equation which follows by demanding O s [n s ] 
to be stationary with respect to H" s : 

A„ f 1 



n - = l T J„ WTKTK - (64 > 

or, equivalently, with m 2 = rh 2 D + II S , 

ml = m|,- -g 2 Tm s . (65) 

This is a second-order algebraic equation for m s , whose r.h.s. is recognized 
as the expansion of the integral in the general gap equation in four dimen- 
sions, eq. (|56|l . up to terms linear in m/T. Note that, to this accuracy, there 
is no ultraviolet divergence at all, so the coupling constant A (or g 2 ) in the 
equations above is the renormalizcd one. 

The free energy (|63|l can be easily evaluated in terms of the solution m s 
to the above gap equation. One obtains [P s = —Cl s /V): 

Note that this contribution starts at order g 3 , but includes terms of all 
higher orders, via m s . The sum P = P( 2 ) + Pa 01 ec l s - and should 
be compared to the fully self-consistent result in eq. (|61|l . Clearly, these two 
approximations coincide up to order g 3 , but deviate in higher orders, with 
the deviations restricted, however, to contributions due to the hard modes. 
In Figs. |H1 and El we also compare the numerical solutions to the full gap 
equation fBl^l and its "soft" approximation and also the corresponding 
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predictions for the pressure. We see that these predictions are extremely 
close to each other, for both m and P, which confirms that the higher order 
effects associated with the hard modes are small indeed. 



2.3.3. The entropy 

We have mentioned already that the stationarity property l|52|) allows one 
to simplify the calculation of the entropy: 

S = -d(n/V)/dT (67) 

or, more generally, of any first partial derivative of the pressure. Indeed, 
the thermodynamic potential l|53|l depends on the temperature through 
the statistical factors n(u>) and the propagator, or equivalently its spectral 
function p. Because of cq. (1521 the temp erature derivative of the spectral 
density cancels out and one obtainjniEol 

f d A k dn(uj) , „ i . 

+ /(£)< dn dT fc) ReD{ ^ k) + s> (68) 



with 



, _ d(T$) f d 4 k dn(oj 



S' = — 



dT 



d J (2tt) 4 dT 



RenimZ?. (69) 



Loosely speaking, the first two terms in eq. (|68|l . which have a one-loop 
structure, represent the entropy of 'independent quasiparticles', while S' 
accounts for a residual interaction among these quasiparticles^' (see also 
the discussio n below) . For 2-loop f skeletons we have the following additional 
simplificatiorpH2S125l 

S' = 0, (70) 

a property which holds independently of the self-consistency condition Q52JI ■ 
The formula l|68(l . together with S' = 0, makes the entropy a privileged 
quantity to study the thermodynamics of ultrarelativistic plasmas. Let us 
discuss some of the attractive features. 



'For the scalar field theory with </> 4 self-interactions, the property 1701 holds also for the 
3— loop skeleton in eq. 1491 . 
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First, we note that, for a finite function H(oj, k), the integrals in eq. (|68|) 
with S' = are manifestly ultraviolet finite, since dn/dT vanishes expo- 
nentially for both to ~ > ±oo. This reduces the problem of the UV renor- 
malization of the 2-loop ^-derivable approximation for the entropy to that 
of the renormalization of the gap equation for the 1-loop self-energy. In 
general, this is still a difficult task (because this is a non-perturbative ap- 
proximation), but as we have seen in the previous subsection, progress 
can be made at least in spec ific cases. In particular, it has been recently 
demonstrated, in Refs.l^SEZI that for scalar field theories at least, the self- 
consistent approximations are renormalizable to all orders, i.e., they can be 
made finite by adding temperature-independent (vacuum) counterterms to 
the Lagrangian, which are themselves determined in a self-consistent way. 

Furthermore, as already mentioned, the entropy has a more direct quasi- 
particle interpretation than the pressure. To see this more explicitly, use the 
following identity: 

/ Im n \ 

Im logD -1 ^,*;) = arctan - — - 7re(w)0(- Re Z?" 1 ), (71) 

y Rc D J 

to rewrite S as S = S po i e + S damp , with 

S pole = I ^^W M-Rea-Vfc)) 

("^31 {C 1 +nfc)log(l + n k ) - n k logn k j, (72) 

where n k = n(sk), with e k solution of Re_D _1 (w = s k ,k) = 0. This is 
the on-shell energy of the dressed single-particle excitations in the system, 
or 'quasiparticles'. The quantity S po i e is the entropy of a system of 'non- 
interacting' quasiparticles, while the quantity 



S, 



damp 



d 4 k dn(ui) 
(2tt) 4 dT 



ImTI Rc D — arctan 



Imn 
ReTF" 



(73) 



which vanishes when Imll vanishes, is a contribution coming from the 
continuum part of the quasiparticle spectral weights. We see that for the 
entropy, unlike for the pressure, the main effects of the interactions come 
via the change in the spectrum of the excitations, and therefore can be 
taken into account by simply dressing the propagator. 

The above formulae, which are generic (i.e., they hold for any field the- 
ory, including QCD), become even simpler when the 2-loop self-consistent 
approximation is applied to the scalar theory with 4 self-interactions. 
Then, the self-energy reduces to a local mass term, as we have seen in the 
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previous subsection, and therefore S = S po i e , with S po i e given by eq. (|72*|) 
with e 2 = k 2 + to 2 and m determined by eq. (|59[1 . That is, in this approx- 
imation, the entropy of the interacting scalar gas is formally identical to 
the entropy of an ideal gas of massive bosons, with mass m. As emphasized 
after eq. (jfilj) . such a simple interpretation does not hold for the pressure. 



2.3.4. Approximately self-consistent calculations of the entropy 

In this subsection, the calculation of the self consistent 2- loop entropy will 
be considered from a new perspective, which will suggest a different strategy 
for implementing the self-consistency in an approximate way. This strategy 
will be especially useful in view of subsequent applications to QCD, where 
a fully self-consistent calculation seems prohibitively difficult, and not even 
desirable (see Sect. 14 .2|) . 

It is first instructive to see how the perturbative effects of order g 2 
and g 3 get built within the self-consistent 2-loop entropy, as given by (cf. 
eq. C2J): 

s{m) --Jw^^bP lm log(fc2 - " 2 + m2) 

2tt 2 T 3 m 2 T m 3 m , , , s 

= + +0(m 4 T). (74 

45 12 12tt V 1 ' y ' 

The first two terms in the expansion in the second line come from inte- 
grating over hard momenta (k ~ T) within the loop integral in the first 
line, while the third term, non-analytic in to 2 , comes from the soft loop 
(k m). 

The expansion in eq. (|74Jl must be considered simultaneously with the 
perturbative expansion of the solution to to the gap equation l)59|l . To 
order g 3 , this is the same as the standard perturbative expansion, eq. J23Jl: 
to 2 ~ g 2 T 2 (\ — Sg/ff). 

If one replaces to in eq. (|74[) by the LO thermal mass rhj) — gT, then the 
resulting approximation for S reproduces the perturbative effect of order 
g 2 , but underestimates the correction of order g 3 by a factor of 4. But to 
order g 3 , the thermal mass itself must be corrected, as discussed above. 
By replacing m 2 — > rh 2 D + 8m 2 D in the second order term of eq. (|T4p , the 
correction 8m 2 D to the thermal mass generates the 3/4 of the plasmon effect 
missing. To conclude, the order-g 3 contribution to the entropy comes up in 
a subtle way: 1/4 of it comes directly from the soft loop in eq. I|74|l . with 
to replaced by the LO screening mass rhp, while the other 3/4 comes from 
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2 2.5 3 

g(fi = 2irT) 

Fig. 10. Comparison of the fully self-consistent 2-loop entropy with the perturbative 
results to order g 2 and g 3 (long and medium dashed lines, resp.). The short-dashed line 
between the order g 2 result and the full line corresponds to the "HTL approximation" 
where the 2-loop entropy is evaluated with just the HTL mass; the short-dashed line 
almost on top of the full line corresponds to the "NLA approximation" using l76l . 



the hard loop, via the NLO correction to the dispersion relation of the hard 
particles. 

Of course, within the present context, where the self-energy is indepen- 
dent of the external momentum, there is no distinction between the screen- 
ing mass and the quasiparticle mass: both are solution to the same gap 
equation (|59|l . But in QCD, these are distinct quantities already at lowest 
order, and this difference will play an important role in the construction of 
approximations in Sect. 01 Thus, although this may look at first somewhat 
artificial here, we shall introduce different notations for the two masses. 
We shall keep the notation mo for the screening (Debye) mass, and denote 
the hard quasiparticle mass as moo', the latter characterizes the excitation 

. Thus 



spectrum at large momenta, typically of the form Ek = \Jk 2 + t 
the mass m 2 in the second order term of eq. I|74|> should be read as m 
while m 3 in the cubic term is m 3 D . That is, we rewrite eq. I|74(l as 



2 

oo ' 



S(m) 



2n 2 T 3 
45 



12 



12tt' 



(75) 



In perturbation theory, mjj and coincide with each other up to 
order g 3 , but differ in order g 4 and higher. 8 Within the 2-loop ^-derivable 
approximation, they coincide if computed exactly, as already mentioned. 



g Note that only the corrections to mo can be calculated by means of dimensional re- 
duction. Dynamical masses in general req uire the resummation of the nonstatic hard 
thermal loops (for simple examples see Ref. 



February 1, 2008 20:44 



WSPC/Trim Size: 9in x 6in for Review Volume 



tdqgp 



Thermodynamics of the High- Temperature Quark-Gluon Plasma 33 

But if the self-consistency is imposed only in the soft sector, the corre- 
sponding approximations may a priori look different. We have seen earlier 
that self-consistency in the soft sector is enough to obtain an excellent ap- 
proximation for the pressure, and we want to verify here the status of such 
an approximation for the entropy. In fact we shall examine the accuracy 
of two non-perturbative approximations to the entropy, that we shall later 
implement in the case of QCD. In both approximations, the integral in the 
first line of eq. H74|) is evaluated exactly (i.e., numerically), but for values 
of m which are only approximate solutions to the gap equation. 

i ) In the 'HTL approximation' we set m = tHd- The ensuing approxi- 
mation for S includes the correct LO effect of the interactions, of order g 2 , 
but only 1/4 of the plasmon effect. In Fig. ^3 the resulting approximation 
is compared to the exact (i.e., fully self-consistent) result, and also to the 
weak coupling expansion. 

ii ) In the 'next -to-leading approximation' (NLA), the mass Woo of the 
hard particles is determined with NLO accuracy, so as to restore the correct 
plasmon effect in S. That is, 

m lo — > m lc ~ (3/v)g 2 Tm D , (76) 

where mr> is the self-consistent mass, given by the solution to the gap 
equation l|65[l . Of course, the mass thus obtained coincides with the 
Debye mass mc, but this 'degeneracy' is specific to the present scalar model, 
as we have already emphasized. 

The quality of the NLA based on eq. I|76|) can be appreciated from its 
comparison, in Fig.^3 with the "exact" (fully self-consistent) result, and 
also with the HTL approximation alluded to before. Although the HTL 
approximation is seen to deviate from the fully self-consistent result at large 
coupling, it remains smaller than the ideal-gas result despite the presence of 
a quarter of the plasmon effect, which, when treated purely perturbatively, 
would result in a nonmonotonic behaviour leading to an excess over the 
ideal-gas result at large coupling. 

The ^-derivable 2-loop approximation in fact becomes "exact" in the 
large- N limit, where N scalars are interacting in an 0(iV)-invariant manner 
according to ((j> 2 ) 2 . In this limit the exact result is ft- independent, while all 
approximations to it are renormalization scale dependent. However, as can 
be seen in Fig. the agreement is almost perfect for NLA when p, = 2nT 
and varies only slightly when varying this scale, which motivates us to 
adopt jl — 2nT as optimal renormalization point. However, in what follows 
we shall always consider variations by a factor of 2 around this central value 
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to have an estimate of the inherent uncertainty of our approximations. 



3. Quark Gluon Plasma 

Before describing how to overcome the difficulties of perturbation theory 
in the thermodynamics of QCD along the lines sketched for the simple 
scalar models of the previous section, we shall briefly review the results of 
perturbation theory obtained so far, which, at least for vanishing chemical 
potential and high temperature, have been recently pushed to the limits of 
calculability. 



3.1. Perturbative results 

In QCD with TV/ flavors of quar ks, the thermodynamic pressure at high 
temperature has been computecpIE21 j n perturbation theory through 
order g 5 : 



P 



45 



T* 



-30 



32 V 4 \ 

AT \ i3/2 



12 ■ 



1 35 

237.2 + 15.977V/ - 0A13N] + — [ 1 



7T 



TV 



'-4- I In 



1 N f 

33 s 
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2^rT/ 



a, . 

7T 
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6 
495 

~Y 



1/2 



-799.2 - 21.967V/ - 1.9267V 



Nf 
6 



In- 



2ttT 



(t) 



5/2 



In a s 



(77) 



where a s — g 2 (ft) / (4?r) . The coefficient of the In a s term, the last in the 
pressure at high T and vanishing chemical potentia l th at is calculable in 
perturbation theory, has recently been determined a 



p 



8ir u 



rp-i 



I 6 In ( 



45 
1485 



1134.8 + 65.897V/ + 7.6537V| 



TV, 



1 - —Nf] In — 
33 1 / 2vrT 
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(78) 
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Here ft, is the renormalization scale parameter of the MS scheme and 
a»(U) is the running coupling whose beta function is known to four- loop 
ordeJ^I though we shall restrict ourselves to its two-loop version in the 
following, because lattice determinations of T c /Aj^g usually employ only 
the latter. 

Evaluating the perturbative expression to order g 2 , g 3 , g 4 , and g 5 gives 
results which do not show any sign of convergence, as depicted in Fig. 1111 
but rather show increasing ambiguities due to the dependence on the renor- 
malization point, signalling a complete loss of predictive power. (When the 
<7 6 In g-teim of Eq. (|7H|l is included with a suitably adjusted constant under 
the log, the latter can at least be chosen such that the lattice results are 
qualitatively reproduced^.) 

P/Po 
1.4 

1.2 

1 

0.8 
0.6 

1.5 2 2.5 3 3.5 4 4.5 5 

T/T c 

Fig. 11. The poor convergence of thermal perturbation theory in pure glue QCD. The 
various grey bands bounded by differently dashed lines show the perturbative results to 
order g 2 , g 3 , g 4 , and g 5 with MS renormalization point p, varied betwee n ttT and 4-7rT. 
The thick dark-grey line shows the lattice results from the Bielefeld groupS5| the lighter 
one behind that of a more rece nt lattice calculation using an RG-improved action from 
the CP-PACS collaboration 

However, the poor convergence behaviour of perturbation theory is not 
very different from that we have observed above in the case of simple scalar 
4> 4 theory. In the following we shall review the attempts to overcome this 
impasse by reorganizing perturbation theory on the one hand by a variant 
of SPT and on the other by an approach which is guided by ^-derivable 
approximations for entropy and quark number density. In both cases, the 
resulting expressions go beyond strict perturbation theory in that higher- 
order terms are kept unexpanded, and the hope is that strictly nonpertur- 
bative contributions that appear in the pressure of hot QCD at order g 6 
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make up only a small, negligible part of what can be captured already by 
an effective quasiparticle description as soon as the temperature is a few 
times the transition temperature. 

When compared with the perturbative results, it is remarkable that 
the next-to-leading result to order g 2 fares rather well at temperatures 
> 2T C , though the higher-order results prove that perturbation theory is 
inconclusive. On the other hand, simple quasiparticle models which describe 
the effective gluonic degrees of freedom by 2N g (N g = N 2 — \) scalar degrees 
of freedom with asymptotic masses taken from a HTL approximation can 
be used quite succe ssfully to model the lattice data by fitting the running 
coupling 



3.2. Lattice results 

The thermodynamic potential has been studied on 4-dimensional lattices 
both for pure-glue QCD and QCD with a number of comparatively light 
(but still massive) quarks by a number of groups. 

In the case of pure-glue QCD, the temperature dependence of the ther- 
mal pressure has been studied in great detail on lattices with varying tem- 
poral extent up to 8 x 32 3 for temperatures up to about i.5T r The con- 
tinuum extrapolation using the standard Wilson action of Ref. i s given 
in Fig. 1111 by the darker-grey line with the width of the line roughly rep- 
resenting the estimated error. Similar studies have been perfo rme d using 
improved lattice actionJSZESI and the results obtained in Ref. ESI U sing a 
renormalization-group improved lattice action are included in Fig. 1111 by 
the lighter-grey line. The latter are systematically higher by about 5 ... 2% 
for temperatures 2 . . . 3.5T C compared to the results of Ref. 

As can be seen in Fig. 1111 the pressure drops sharply as the transition 
temperature is approached from above. Not shown in the figure is that 
the pressure in fact drops to virtually zero at and below the transition 
temperature. There the pressure is exponentially suppressed by the rather 
large masses of glueballs. 

In the presence of (light) quarks, lattice calculations are much more 
difficult to perform and the continuum extrapolation is correspondingly 
more involved and uncertain. However, through the use of i mproved lattice 
action, considerable progress has been made recently ^^^^ l When rescaled 
to the ideal-gas values, the results, as a function of T/T c with the respective 
different values of T c , the pattern of the pressure of hot QCD with quarks 
is quite similar to that of pure-glue QCD, and there is moreover a rather 
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weak dependence on the number of flavours if all of them arc light. 
3.3. Dimensional reduction 

Dimensional reduction in hot QCD (at zero che mical potential) leads to an 
effective three-dimensional Lag 

L E = \teF? j + tv[D h A f + m 2 E tr A 2 + ^\ E (tr A 2 ) 2 + ... (79) 

where the parameters are determined perturbatively by m atch in j^^U I In 
lowest order one has a dimensionful coupling g 2 E — g 2 T an 



m| = (1 + N f /6)g 2 T 2 , X E = ^rV T ' ( so ) 



though Xe starts to contribute to the pressure only at order g 6 . At this 
order, however, the self-interactions of the massless magnetostatic gluons 
start to contribute, but these contributions are inherently non-perturbative 
because the three-dimensional theory for the zero modes A (x) is a confining 

theory !n™ 

The thermal pressure of the 4-dimensional theory can be decomposed 
into contributions from the hard modes ~ T, calculable by standard per- 
turbation theory, and soft contributions governed by (|79() which involves 
both perturbatively calculable contributions up to order g 5 T 4 and the non- 
perturbati ve o nes starting at order g 6 T 4 . 

In Ref. ESI the effective theory based on (|79|) has been used to organize 
and r eprodu ce the perturbative calculation of the thermal pressure (I77|l of 
Refs. Ill This turns out to be particularly elegant when dimensional 
regularization is used to provide both the UV and IR cutoffs of the original 
and effective field theories. 

T o order <? 4 , the contribution of the hard modes can then be written 
aJSl 

8tt 2 _ 4 U. . 21 , T \ , -15 A . 5„\a, 



P haId = -T^n + -N f + —n + -N f 



244.9 + 17.247V/ - 0.4157V? + 135 1 + —± In 



6 J 2ttT 



'"^'X 1 - 5"') < 81 > 



In the first logarithm the dimensional regularization scale p, is associated 
with regularization in the infrared and thus has to match a similar loga- 
rithm in the effective theory, whereas the second logarithm is from UV and 
involves the first coefficient of the beta function. 
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Indeed, calculating the pressure contribution of the soft sector described 
by 1)79(1 in dimensional regularization gives, to three-loop order (neglecting 
Xe -contributions ) 



P soft /T = —ml - -At ( 4 In + 3 ) glml 
9/89 11, „ 1 



8.3 v 24 6 ™+r* ( 82 ) 

All the contributions to the pressure involving odd powers of g in 1(77(1 (as 
well as part of those involving even powers) are coming from the soft sector. 
Inserting the leading-order value l(80() for itie gives the QCD pressure up 
to and including order g^lng; to obtain all the terms to order g 5 , next-to- 
leading order corrections to the m ^- par ameter have to be obtained by a 
matching calculation as given in Ref. ™l. 

In order to obtain the g 6 In ^-contribution 1)78(1 one also needs g E to 
order g 4 (given in Ref. I^J) and above all the four-loop contribution of the 
effective theory 1(791) which has recently been calculated analytically a 



p (4) /T _ N (Nglf 



43 157tt 2 1 



12 768 



In 

9e 



43 491tH 



4 768 



ln-£=- 

m E 



(83) 



up to a constant c that is strictly nonpcrturbative and needs to be de- 
termined by three-dimensional lattice calculations. Such calculations have 
[T7I 

been carried out recently^, but they depend on an as yet undetermined 4- 
loop matching coefficient. At the moment the conclusion is that it is at least 
not excluded that the lattice results based on dimensional reduction can be 
matched to the full four-dimensional results at temperatures of a few times 
the transition temperature. For this reason, the most reliable results on the 
thermodynamics of hot QCD (particularly for pure-glue h QCD) remain to 
date the four-dimensional lattice data. 

An interesting feature of the perturbative calculation when organized 
through the dimensionally reduced effective theory ((791) is that the large 
scale dependences of strict perturbation theory can be reduced when the 
effective parameters are not subsequently expanded out when their per- 
turbative results are inser ted in dH2J 63 7 ^ What is more, considering suc- 
cessively the full two-loop and full three-loop contributions from the soft 



'inclusion of fermions is particularly easy in the dimensional reduction method, which 
makes a full three-dimensional prediction clearly most desirable. 
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sector, the results do not exceed the ideal-gas result (as the strictly pertur- 
bative results to order g 3 and g A do), a nd their scale dependence diminishes 
by going from two to three-loop order!^ 

At three-loop order, it is in fact possible to eliminate the scale depen- 
dence altogether by a principle of minimal sensitivity. The result in fact 
agrees remarkably well with the 4-d lattice results down to ~ 2.5T C . (This 
is also true when the four-loop logarithms are included provided a suitable 
value for the unknown constant c in l)83|) is chosenP^*^) 

While this goes only minimally beyond a strictly perturbative treat- 
ment, it strongly suggests that perturbative QCD at high temperature, 
when supplemented by appropriate resummation of soft physics, is capable 
of quantitative predictions at temperatures of interest. 



4. Thermodynamics of HTL quasiparticles in QCD 
4.1. HTL- screened ■perturbation theory 

In the simple example of scalar <p A theory, we have seen that a promising 
strategy for improving the poor convergence of perturbation theory may be 
screened perturbation theory based on the reorganization of the Lagrangian 
according to (|37|l . The corresponding resummation changes the UV struc- 
ture of theory at any finite order of perturbation theory, but removing 
the additional UV divergences by some minimal subtraction does seem to 
lead to results which avoid the bad behaviour of conventional perturbation 
theory. 

This method can be extended to QCD by employing the gauge invari- 
ant HTL effective actior P I ^ I 1 ^ instead of a simple mass term. The HTL 
effective action is nonlocal and involves nontrivial vertices in addition to 
self-energies, but it is uniquely determined by only one ma ss parameter 
for gluons (and a further one for quarks, if any-jZHZBESEB This variant 
of screened perturbation theory has been termed HTL perturbation the- 
ory (HTLPT) 1 by their authors, but it should be noticed that it is meant 
to transcend perturbative HTL resummation, where hard propagators and 
vertices with hard propagators attached remain undresseo . 

In HTLPT the QCD Lagrangian is improved by writing 

£ncn = £ 



QCD = •'-QCD 



(1 - J)£ H tl (84) 

jrv t 

recall that this method transcends the original HTL resummed perturbation theory*!. 
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with the standard HTL Lag 

rangiarPEZEl 

£htl = ~\ m ° trF - (www) F ^ + m2 ^ (dfey ) ( 85 ) 

where = (1, v) is a light-like four-vector and (. . .) represents the average 
over the directions v. Since this only serves as gauge invariant generalization 
of the thermal mass term for SPT, i.e. a variationally improved perturbation 
theory, the parameters niD and M need no longer be restricted to their 
leading-order HTL values (given in eqs. I|91|l and (|94|l below). Technically, 
HTLPT corresponds to using S as an expansion parameter that is set equal 
to 1 event ually . 

In Ref. HTLPT has been used to calculate the pressure of QCD at 
one-loop level. This corresponds to calculating the thermal pressure of a 
modified free Lagrangian l|84|) with S = 0. If (but only if) the mass param- 
eters mo and M are identified with their leading-order HTL values (|91|) 
and l|94|) , this includes the plasmon term ~ g 3 of the pressure correctly, but 
over -inc ludes the leading-order interaction term ~ g 2 by a factor of 2. (In 
Ref. 121 the over-inclusion factor was in fact larger because of an inconsis- 
tent usage of dimensional regulariza tiorl^l Corrected numerical results at 
one-loop level are contained in Refs. E2E0I ) 

Starting at two-loop order, it becomes possible to determine mo and M 
by a variational principle. The two-loop pressure in HTLPT is in fact almost 
prohibitively complicated because of the necessity to include th e non- local 
HTL vertices along with the nontrivial HTL propagators. Refs. ISHEDl have 
simplified this task by expanding the resulting sum-integrals in powers of 
mo/T and M/T such as to include all terms through order g 5 if tud/T and 
M/T are taken to be of order g. (These mass expansions have been checked 
to be fairly accurate in the applications of intereslr^.) When determining 
mo/T and M/T through variational equations, the two-loop HTLPT result 
is equivalent with the perturbative result (|77|l through order g A log(g), with 
deviations setting in at order g 4 . The thermal mass parameter rajj coincides 
with the perturbative result for the Debye mass at leading order, though not 
beyond, which is not too surprising since it is just a variational parameter 
and not directly related to the full electrostatic propagator; the quark mass 
parameter M when determined by a variational principle comes out as ~ g 
at weak coupling, however with a coefficient that differs from leading-order 
perturbation theory and is in fact subtraction-scale dependenlP"2 

When evaluated numerically, the HTLPT results, like those of SPT in 
the scalar case, turn out to be a clear improvement over the conventional 
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0.95 : 



0.85 




Fig. 12. The two-loop DRSPT resultf^lfor pure- glue QCD (gray band) in comparison 
with the two-loop HTLPT result of Ref. EH (dashed lines) and the "NLA" result of 
Ref. 1^1 described in Sect. H~4l (dash-dotted lines), all with fl varied around 27rT by a 
factor of 2. 



perturbative results in that (at least for all T > 1.5T C ) both the one-loop 
and the two-loop results are lower than the ideal-gas result while including 
the plasmon effect ~ g 3 which in conventional perturbation theory (see 
Fig. leads to results much above the ideal-gas value in contradiction 
with all lattice results. Moreover, the one-loop and the two-loop results 
do not deviate too strongly from each other so that they seem to have 
satisfactory convergence. On the other hand, they fall short of describing 
lattice results quantitatively but remain too close to the ideal-gas value and 
account for only about half of the deviation of the lattice results from the 
latter. 

A hint on the possible reason for this apparent deficiency of HTLPT 
in the application to the thermodynamics of QCD comes from a simpler 
variant of SPT, namely when the latter is implemented only in the dimen- 
sionally reduced sector (j79(l by rewriting l|79|l trivially as 



1 



1 



1 



1 



£eqcd = i ^^ + -(A^o) Q (AAo) a +-(m^+W)AgA a ~-WA^g.(86) 

Here m 2 E is determined by perturbative matching, whereas the additional 
5m 2 is to be fixed by a gap equation (i.e., a principle of minimal sensitiv- 
ity). Like any SPT, this dimensionally reduced SPT (DRSPT) introduces 
additional UV singularities at finite loop order, which should be minimally 
subtracted to allow comparison with HTLPT. 

Like HTLPT, this formalism is gauge invariant, but very much simpler 
and it avoids the intermediate modification of hard modes (where ordinary 
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perturbation should be accurate anyway). In Ref. ^ it has been recently 
shown that the two-loop result of DRSPT is much closer to the lattice 
results (reproduced in Fig. I12|) . Because fermions are integrated out in 
DRSPT, this also avoids the problem HTLPT appears to have with a gap 
equation for the fermionic mass. (On the other hand and in contrast to 
HTLPT, DRSPT is clearly restricted to static quantities.) 

Unfortunately, at three-loop order DRSPT gives rise to a gap equation 
for Sm 2 which does not connect to perturbation theory any lo nge r so that 
it cannot be viewed as an improvement of perturbation theory!^ 

4.2. 2PI formalism in gauge theories 

In the following we shall describe a different approach in which the 2PI 
formalism in terms of full propagators is used to construct approximations 
which avoid the UV problems and ambiguities of SPT and HTLPT. This 
turns out to amount, at the two-loop level, to a description of the thermo- 
dynamics of hot QCD in terms of weakly interacting quasiparticles, whose 
nontrivial spectral data serve to incorporate (at least in leading and sub- 
leading orders) the strong interactions of the elementary degrees of freedom. 

In QCD, the thermodynamic potential is a functional of the full gluon 
(D), quark (S), and Faddeev- Popov ghost (D g h) propagators, 

PSI[D, S, D gh ] = 1 Tr log IT 1 - Tr logS" 1 - Tr log^ 1 

-~TrILD+ TrVS+ TrIL gh D gh + 3>[D,S,D gh },(87) 

where Tr now includes traces over color indices, and also over Lorentz and 
spinor indices when applicable. The self-energies for gluons, quarks and 
ghosts are denoted by n, E and U g h , respectively. 




(d) 



Fig. 13. Diagrams for $ at 2-loop order in QCD. Wiggly, plain, and dotted lines refer 
respectively to gluons, quarks, and ghosts. 



In the following we shall consider a two-loop ^-derivable approximation 
which amounts to keeping the diagrams displayed in Fig. ED This again 
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leads to the important simplification that the entropy is given by a com- 
paratively simple expression 



with vanishing interaction term 5' (cf. (|69(l ). and this also holds true for 
fcrmion number densities: 



where the completely analogous interaction term M' is equally absent at 
the two-loop level. 

While ^-derivable approximations are symmetry preserving as concerns 
global symmetried^, local gauge invariance is not respected, so this approx- 
imation introduces gauge dependences as long as vertices are not dressed 
in a self-consistent manner in line with the self-energies J . These gauge de- 
pendences are actually parametrically suppressed at the stationary poinlP^J 
where they enter only at twice the order of the truncation. However, in what 
follows we shall construct further approximations based on the gauge invari- 
ant hard thermal loops, and this will entail manifest gauge independence 
in our applications. In these approximations, where the self-consistency of 
the two-loop ^-derivable scheme is only approximately realized, but such 
that the deviations from full self-consistency are of the order of three-loop 
corrections, only the two physical polarizations of the gluons contribute 
and the Faddeev-Popov ghosts just cancel the unphysical gauge degrees of 
freedom of the former so that they can be simply dropped. 

4.3. Approximately self- consistent entropy and quark 



As discussed above, conventional thermal perturbation theory exhibits very 
poor apparent convergence as soon as collective phenomena such as the 
plasmon effect are included, which involve odd powers in the coupling. The 
self-consistent entropy and density functionals at two-loop order also con- 
tain the plasmon effect, but together with higher-order terms that ordinarily 
would have been discarded by a strictly perturbative expansion in g because 




(88) 




(89) 



density 



j(Dn a formal level this has been worked out in Ref. IHOJ 
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otherwise the renormalization programme could not have been carried out. 
By contrast, the two-loop entropy and density functionals 1(88(1 and 1)89(1 are 
UV finite expressions and can be evaluated in a nonperturbative manner 
using the high-temperature/density approximation of the gluon and quark 
propagators, which are finite at leading and next-to-leading order in HTL 
perturbation theory. 



4.3.1. HTL approximation 

As a first approximation we consider the two-loop entropy and density 
functionals evaluated completely using HTL (HDL) propagators. From the 
discussion of the scalar model we anticipate that this includes correctly the 
leading-order interaction contribution cx g 2 in the entropy (in contrast to 
the pressure) . k 

But unlike the scalar 4 model, where the HTL self-energy is a simple 
constant mass, cf. 1(18(1 . the HTL self-energies for gluons and quarks contain 
two non-local structure functions each. The gluon HTL self-energy contains 
transverse and longitudinal components given by 



where 



fl L (w, k) = m 2 D 



m 2 D = (2N + N f ) 



u 1 

— log 

2k B 


UJ 


+ k 




- k 


uj 2 - 


k 2 


n L 


k 2 





g 2 T 2 



N, 



2tt 2 



(90) 



(91) 



is the (leading-order) Debye screening mass. 1 As is well knowJSEl the 
gluon propagator dressed by these self-energies has quasiparticle poles at 
u>T,L{k) with momentum-dependent effective masses and Landau damping 
cuts for \uj\ < k. When k ^> mn, the pole corresponding to the collec- 
tive longitudinal excitation has exponentially vanishing residud^, whereas 



k A different procedure for including the full HTL propagators has been proposed by 
PeshieJH2| who considers a (somewhat ad-hoc) modification of the <I> functional such 
as to have correct leading-order contributions and UV finitcness when full propagators 
are approximated by HTL ones. This turns out to be equivalent to our results when 
expanded out to order g 3 . 

'For simplicity we shall write out our formulae for only one value of the chemical poten- 
tial; the generalization to several different chemical potentials is straightforward. 
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transverse excitations tend to uj — > \Jk 2 + to 2 ^, where the asymptotic mass 
is given by 

m 2 00 =fl T (k,k) = ^m 2 D . (92) 

The (massless) quark propagator at finite temperature or density is 
split into two separate branches of opposite ratio of chirality over helicity 
with propagators A± = [—uj + k ± X±] _1 . In the HTL approximation, the 
respective self-energies readEHl 

p . . M 2 ( uj T k , uj + k\ . , 



where (C/ = (TV 2 - l)/2iV) 



7T 



M = IT 2 + ^ . (94) 



The dressed fermion propagators have quasiparticle poles at w± (fc) and 
Landau damping cuts. At large k and positive frequency uj + — > ^Jk^ + M 2 ^ 
with asymptotic mass = 2M 2 ; the additional branch cj_ which has 
no an alogue at zero temperature and density has exponentially vanishing 
rcsidui 



Although the various HTL self-energies constitute a leading-order result 
only under the condition of soft momenta and frequencies, the result for the 
asymptotic masses applies also for hard moment d^ l ^ l 

Using these expressions in the entropy formula l|88l) this takes care of all 
contributions of order <? 2 , but only part (< 1/4) of the plasmon term ~ g 3 . 
However, it also contains infinitely many higher-order terms which despite 
being incomplete may help to get rid of the pathological behaviour of the 
perturbation series truncated at low orders in g. Among such higher-order 
contributions is for instance a g A contribution reading (for pure glue) 



4 4 k = -N g7 z3z ( log ^- + 1-55 . . . ) (95) 



T 

16n' 2 T \ ° rhry 

involving a g 4 log(l/g)-term, whose coefficient in pure-glue QCD is 1/12 
of that of the complete perturbative result. 111 By co ntrast, simple mas- 
sive quasiparticle models such as those used in Refs. lalLagl d no t have 
a g 4 log(l/<7)-term in the entropy at all. 



m The correct coefficient will be restored by 0(g 4 log(l/g)T 2 ) corrections to m^, whereas 
the constant under the logarithm also receives three-loop contributions. 
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Shtl/Ssb 
1 

0.9 
0.8 
0.7 
0.6 



0.5 1 1.5 2 2.5 3 

rhn/T 

Fig. 14. HTL entropy per gluonic degree of freedom normalized to its Stefan-Boltzmann 
value as a function of the Debye mass mo(T, /i)/T. The full line gives the com- 
plete numerical result; the dashed line corresponds to the perturbative result to order 
(m D /T) 3 ~ g 3 . The dotted line gives the entropy for scalar degrees of freedom with 
momentum-independent mass m = rriao = m,£, / v2; its perturbative approximant is 
given by the dash-dotted line. 



In Fig. 1141 the numerical result for Shtl/Ssb in the case of pure glue is 
given as a function of rho/T, which is the only independent parameter. The 
HTL result is depicted by the full line and is found to be a monotonically 
decreasing function of tjid/T. If this were expanded in powers of rho/T 
and truncated beyond (mu/T) 3 ~ g 3 (dashed line in Fig. I14|). this property 
would have been lost at m£>/(27rT) w 1/3, where one might still expect a 
sufficiently clear separation of hard and soft scales which is a prerequisite 
of the HTL approximation. 

The numerical result for Shtl is in fact very close to a simple mas- 
sive quasiparticle model with entropy 2N g SsB('nioo), represented by the 
dotted line in Fig. This is due to large cancellations between the lon- 
gitudinal and transverse contributions to the HTL entropy. However, quite 
remarkably one has Shtl < 2N g SsB(moo) even though the latter contains 
30% less of the plasmon term ~ g 3 , which would be expected to work in 
the opposite direction. This is further emphasized by comparison with its 
perturbative approximation given by the dash-dotted line. 

In Fig.EDthe HTL entropy of 1 quark degree of freedom at zero chemical 
potential is displayed as a function of the fermionic plasma frequency M. 
Like in the pure-glue case, this contains also an (incomplete) g 4 log(l/<7) 
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contributions that is not present in simpler quasiparticle models: 

S [ f %TL = NN f^f ( lQ g § + °- 22 • • •) (96) 

Sf,HTL does not contain any contributions to the plasmon term ~ g 3 , which 
entirely come from NLO corrections to . There is therefore rather little 
deviation of the full numerical result for the fermionic contribution (solid 
line) from the perturbative one truncated above order g 2 (dashed line). 

Compared to the entropy of a simple massive fermionic quasiparticle 
model 5} 0) (Moo) (dotted line in Fig. HSJ), one finds extremely good nu- 
merical agreement, which again takes place only after adding up all the 
quasiparticle ((+) and (— )) and Landau-damping contributions, however. 

Sf.HTL/Sf* 1 
1 

0.9 

0.8 
0.7 
0.6 



0.2 0.4 0.6 0.8 1 1.2 1.4 

M/T 

Fig. 15. HTL entropy per quark degree of freedom at fj, = normalized to its free 
value (solid line), the corresponding perturbative order-g 2 result (dashed line), and the 
entropy of a quark with constant mass Moo (dotted line). 

In the case of zero temperature but high chemical potential, the quantity 
of interest is the quark density. This does not contain any plasmon term 
~ g , but rather ~ <? 4 log(l/g). The HDL approximation does contain some 
though not all of this term: 

J^hL = NN f^ ( lQ g f { + °- 35 • • ■) (97) 

Order-<? 2 log g contributions to the asymptotic masses of the quark and 
gluon quasiparticles, still within the framework of the 2-loop ^-derivable 
approximation, are responsible for the remaining contribution to the coef- 
ficient of the g 4 log g-term, while the coefficient under the logarithm also 
receives 3-loop contributions. 
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In Fig.^|the numerical result for Mrdl is given as a function of M//i 
(full line) and compared to that of a simple quasiparticle model 



= — ^ (/i 2 — M^)3#(/i — Moo) 



(98) 



as well as to a perturbative approximation truncated beyond [M / ' \i) 2 ~ g 2 . 
Remarkably, the full HDL result drops to zero at almost the same value 
as the simple quasiparticle model. However, the former becomes nega- 
tive thereafter, showing that the HDL approximation is breaking down at 
-^oo/V > 1 at the latest. 



M„ DL /M^ 

1 




0.6 0.7 

Fig. 16. HDL quark density at T = normalized to its free value (solid line), the 
corresponding perturbative order-g 2 result (dashed line), and the free quark density of 
a quark with constant mass Moo (dotted line). 



4.3.2. Next-to-leading approximations 

The plasmon term ~ g 3 at high temperatures T 3> Tod becomes com- 
plete only after inclusion of the next-to-lcading correction to the asymp- 
totic thermal masses and . These are determined in standard HTL 
perturbation theory through 



= Re5U T {uj = k) 



Re( ~\~ ^aa/IaaaaaIaa/ 1 ~\~ jv\ivw\w — A 



(99) 
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where thick dashed and wiggly lines with a blob represent HTL propagators 
for longitudinal and transverse polarizations, respectively. Similarly, 

^SMl(k) = SS+(w = k) 

= Re( C X )\ u = k . 

The explicit proof that these contributions indeed restore the correct 
plasmon term is given in Ref. "SI 

These corrections to the asymptotic thermal masses are, in contrast to 
the latter, nontrivial functions of the momentum, which can be evaluated 
only numerically. However, as far as the generation of the plasmon term is 
concerned, these functions contribute in the averaged form 

- 2 _ Jdkk n' BE (k) Re 6IL T (u = k) 

°°~ Jdkkn' BE (k) (W1) 



and similarly 



/ dkkn' FD (k) Re2fc<5£ + (w = k) 



5M L = T, , } ns L ■ ( 102 ) 

00 Jdkkn' FD (k) v ; 

which can be determined analytically , leading to the remarkably simple 
expressions 

dm^ = -l-g^NTmn, (103) 
SMl = --Lg2c f Tm D ,C f = N g /(2N), (104) 

Z7T 

These corrections only pertain to the hard excitations; in contrast to 
the scalar case, in gauge theories the one-loop corrections to the vari- 
ous thermal masses of soft excitations are known to differ substantially 
from 1103(1 or 1(104(1 . For instance, the gluonic plasma frequency at k = 
readJSH Jm^; /m^ sa — Q.18y/Ng, which is only about a third of Sm^/rrr^; 
the NLO correction to the nonabelian Debye mass on the other hand is 
even positive for small coupling and moreover logarithmically enhanced^ 
5m 2 D /m 2 D = +s/W/(2n) x ff log(l/ ff ). 

For an estimate of the effects of a proper incorporation of the next-to- 
leading order corrections we have therefore proposed to include the latter 
only for hard excitations and to define our next-to-leading approximation 
(for gluons) through 



Snla — Shtl 



soft 



(105) 

hard 
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where m?^ includes (|1(J3|) . To separate soft (k ~ mo) and hard (k ~ 2nT) 
momentum scales, we introduce the intermediate scale A = ^^ixTrtiuCK 
and consider a variation of ca = \ ■ ■ ■ 2 as part of our theoretical uncer- 
tainty. 

Another crucial issue concerns the definition of the corrected asymptotic 
mass moo. For the range of coupling constants of interest (g > 1), the 
correction IJm^J is greater than the LO value m^, leading to tachyonic 
masses if included in a strictly perturbative manner. 

In the simple scalar model of Sect. 12.3121 we have seen that the require- 
ment of self-consistency for the soft modes leads to a quadratic gap equation 
for the asymptotic thermal mass, eq. i|76|) . 

In QCD, where the (non-local) gap equations are much too complicated 
to be attacked directly, we simply consider perturbatively equivalent ex- 
pressions for the corrected TOoo which are monotonic functions in g. Besides 
the solution to a quadratic equation analogous to 1)76(1 . namely 

_ 2 g*(N + N f /2)T* g 2 NT_ 

m x = "717™- ( 106 ) 

we have tried the simplest Pade approximant (degree [2,1]), which also 
gives a greatly improved approximation to the solution of the gap equation 
of the scalar (large N) model. In QCD, our final results do not depend too 
much o n w hether we use the Pade approximanJ23^1 or a quadratic gap 
equatiorPS. 

The main uncertainty rather comes from the choice of the renormaliza- 
tion scale which determines the magnitude of the running strong coupling 
constant (chosen to satisfy the renormalization group equation to 2-loop 
order in the following). 

In Fig. El the numerical results for the HTL entropy and the NLA one 
are given as a function of T/T c with T c chosen as T c = 1.14Ajjg. The full 
lines show the range of results for Shtl when the renormalization scale ft 
is varied from nT to AnT; the dash-dotted lines mark the corresponding 
results for Snla with the additional var iati on of ca from 1/2 to 2. The 
dar k-gr ay band are lattice data from Ref. ^° (the more recent results from 
Ref.'^are consistent with the former within error bars and centered around 
the upper boundary of the gray band for T w 3T C and somewhat flatter 
around 2T C ). Evidently, there is very good agreement for T > 2.bT c . 

When Nf > massless quarks are included, the next-to-leading correc- 
tions to the asymptotic fermionic masses ((104JI have the same problem as the 
bosonic ones that for large coupling the asymptotic masses may turn tachy- 
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Fig- 1 7- Comparison of the lattice data for the entropy of pure-glue SU(3) gauge theory 
of Ref. I^^l (gray band) with the range of Shtl (solid lines) and Snla (dash-dotted lines) 
for p, = nT . . . 4ttT and c A = 1/2 . . . 2. 



omc. In Ref. 1^3 they have therefore been treated like the bosonic asymptotic 
masses, by postulating the simplest quadratic gap equation which is per- 
turbatively equivalent to l|ll)4|) . However, there is an obvious asymmetry 
in that only gluonic quantities appear on the right-hand side of (|1U3[I and 
(|104|l . Indeed, treating them in the same manner is not con sistent in the 
limit of Nf — > oo, which has been solved recently in Refs. 1 ^ 1 . To be consis- 
tent with this limit, the fermionic gap equation has to remain linear as it is 
in perturbation theory. We therefore restrict our requirement of a perturba- 
tive equivalent gap equation with monotonic behaviour only to the bosonic 
(gluonic) sector. For the fermions we proceed by treating the rig ht-h and 
sides of eqs. (|103ll and (|104|l in the same manner, which amounts tcE5 

g 2 N f C f T* 



Mi 



1 



■g 2 CfTm c 



(107) 



4 V2vr" 

There is then no negative feedback from the fermion mass itself, it only in- 
herits higher-order terms from the solution to . This does not avoid the 
appearance of tachyonic masses when both Nf and g are large, but it turns 
out that for Nf < 3 one still has a monotonic behaviour of M^ as a function 
of g. In fact, in the case Nf = 3 the solutions to the correspondingly cor- 
rected gap equations happen to coincide precisely with those of the (struc- 
turally rather different) decoupled quadratic eq uati ons adopted in Ref. ^3. 
As a result, the changes on the results of Ref. ^3 are rather minimal for 
phenomenologically interesting values of Nf. (In the limit Nf — > oo, there 
are more important effects, such that the picture of stable quasi-particles 
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Fig. 18. Comparison of the HTL entropy (solid lines) and the NLA estimate (dash- 
dotted lines) for Nt = 0, 2, 3 with the estimated continuum extrapolation of Nf = 2 
lattice data of Ref. (The gray dash-dotted line corresponds to using a quadratic 
fermionic gap equation, which only has an effect for NLA Nf = 2, which is minimally 
higher t han the result obtained with a linear one required for consistency with the large- 
N f limit 93 ). 



Nf = 0,2,3 (NLA) 




N f = 0,3, 2 (HTL) 



seems to have to be abandoned above some large coupling, but up to this 
point the quantitative agreement with the exact result is quite goo 

In Fig. El the (updated) results for the en trop y with Nf < 3 mass- 
less quarks are compared with a recent estimate^ of the continuum limit 
of lattice results for Nt = 2 (gray band), but now with Srtl and Snla 
evaluated for the central choice of p, = 2ttT and ca = 1 (with unchanged 
T c /Aj^-g). When Nf is increased, there are competing effects of larger ther- 
mal masses versus slower running of a s , which result into a rather weak 
dependence of our results on Nf as a function of T/A^g as it is in Fig. El 
This seems to be consistent with the very small Nf dependence observed 
on the lattice (when the thermodynamic potentials are normalized by their 
ideal-gas limit and expressed as a function of T/T c ). 

It should be understood, however, that the simple quadratic gap equa- 
tion l|106(l and its extension to the fermionic sector, eq. (|1()7J) . are at present 
only an estimate of how the next-to- leading order corrections (|99(l and 1)100(1 
in actuality modify the spectral properties of the hard quasiparticles. More- 
over, the former are included so far in the averaged form 1(101[1 and (|102|l . 
A possible refinement is to include the non-trivial momentum dependence 
of l|99[l and l(100[l which we expect to render a separation scale between 
hard and soft quasiparticles superfluous because the corrections and 
(|100(1 appear to grow with momentum so that their incompleteness at soft 
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momenta may be negligible. When doing so we also hope to be able to come 
to a better understanding of the mechanism which provides the nonlinear 
feedback preventing tachyonic masses from appearing. 

4.4. Pressure 

In the above, we have concentrated on the entropy (and quark number 
density) as the thermodynamic potential which is most suitable for a (HTL) 
quasiparticle description. 

Were it not for the trace anomaljESI (Tp = £ — 3P ^ the correspond- 
ing pressure would be simply given by P = (TS + /jM) /4 (in the case of 
massless quarks) . The correct relation is instead provided by an integration 
such as 



Up to a single integration constant, entropy and density evidently con- 
tain all the relevant information. This integration constant, however, rep- 
resents strictly nonperturbative information, because it amounts to fixing 
the bag constant. Since QCD is asymptotically free, but confining at low 
temperature, this constant cannot be fixed by requiring that the thermal 
pressure as obtained from entropy and density vanishes at zero tempera- 
ture. It is also not possible to determine P{T{) from the requirement that 
bothp(a s ) = P(T)/P (T) and s{a s ) = S(T)/S (T) approach 1 in the limit 
of a s — > 0, because the differential equation p(a s ) + \(i{a s )p' '(a s ) = s(a s ) 
with P(a s ) = — Poa^ — 0x<y.g — . . . has as homogeneous solutions 



and this vanishes at a s — together with all its derivatives. 

Instead we can fix the integration constant by comparison with lattice 
data for the pressure at fj, = 0. A simple possibility for our purposes is to 
take Ti = T c and P(T C , 0) = 0. For T substantially larger than T c this does 
not influence much the results derived from entropy and density, since in 
the ratio P/Pq the (bag) constant term decays like T -4 . 

The result obtained from integrating the entropy in the case of pure-glue 
QCD with boundary condition P(T C ) = is shown in Fig.ll2land compared 
with available lattice data as well as the results obtained recently in 2-loop 
HTLPT by Andersen et a/P and in a dimensionally reduced variant of 




(108) 




SPlP. 
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4.4.1. Finite chemical potential 

For finite chemical potential, the entropy and quark density expressions 
satisfy the Maxwell relations 



dS 



dN 
~df~ 



(109) 



only to the same order in g to which they are accurate when expanded out 
perturbatively. 

In simple quasiparticle modelJ^EH one therefore has to allow for a bag 
function (instead of a bag constant) to restore thermodynamic consistency. 
This can in fact be utilized for matching lattice results and mapping them 
in a phenomenological manner to nonvanishing chemical poten tial^ , where 
recently significant progress has been made on the lattice side^l This in- 
cludes by now concrete results on the thermodynamical potentialJ^i which 
have been found to agree well with the simple quasiparticle models^ 21 . 

A similar programme has been carried through with the HTL expres- 
sions for entropy and density in place of those of the simpler quasiparticle 
models with constant m asse s, and the results are, like those of the sim- 
pler models, encouraging^! An extension which includes next-to-leading 
corrections to the asymptotic thermal masses is work in progress. 



4.5. Quark number susceptibilities 

Quark number susceptibilities (QNS) are defined as the second derivatives 
of the pressure with respect to the chemical potentials of the different quark 
flavors: 

dMi d 2 p 

where i, j are flavor indices. 

These quantities are of direct interest for studying the transition to the 
quark-gluon plasma phase as they can be related to fluctuations in con- 
served charges which can discriminate against a purely hadronic phase^^l 

O n the theoretical s ide, there are new results on QNS from lattice gauge 
theory^ ^P^P^P^ , extending them to higher T/T c which can be com- 
pared against conventional perturbative results and those resulting from 
resummation schemes aiming at overcoming the poor convergence of the 
former. 
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4.5.1. Diagonal quark number susceptibilities 

The diagonal QNS, which for degenerate quark flavors are all equal and 
denoted simply by x m what follows, are found to increase sharply at the 
deconfinement phase transition toward a large percentage of the ideal gas 
value xo = NT 2 / 3 (for SU(iV) and massless quarks). 

In perturbation theory, the series expansion of x/xo is now known^S 
through order g 6 ln(l/g) and the coefficient under the log at order g 6 is in 
fact computable in perturbation theory, in contrast to the one appearing in 
the thermal pressure, though not yet available. The first few terms of this 
expansion read 1 Y 

X_, 1 3 N g ( g^ 3 N IN N f ( g^ 



Xo 2N 

>s fi.y log- + 0(g^ (111) 



3N a rnsA 



(for the more leng thy results on the subsequent coefficients we refer the 
reader to Ref.EES). 

In these strictly perturbative results the plasmon term overcompensates 
the leading-order interaction term, much like in the case of the pressure, 
though less severely so (for Nf — 2 the diagonal QNS fall below the ideal- 
gas value only when T > 40T C , and only for T > 700T C does one find a 
growth of x/xo with temperature as observed on the lattice — in the case 
of the pressure these temperatures below which perturbation theory clearly 
fails are orders of magnitudes higher). 

The ^-derivable 2-loop expression (|89() captures both the leading-order 
term oc g 2 and the plasmon term oc g 3 accurately, together with a subset 
of higher-order terms that are usually expanded out in 2-loop perturbation 
theory. 

As additional approximations we may start from the HTL quark density 
Mhtl where Q89JI is evaluated with HTL propagators. Differentiation with 
respect to the chemical potential yields the corresponding approximation 
for the diagonal 11 QNS. 

While this approximation contains the leading-order term oc g 2 accu- 
rately, it does not include anything of the plasmon term. In the QNS, the 
plasmon term arises completely from corrections to the asymptotic thermal 
masses of quarks according to ()10()(l . 



"In the HTL approximation the QNS are purely diagonal. 
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Fig. 19. Comparison of our results for x/xo m the formal limit Nt=0 with two recent 
continuum-extrapolated lattice results for qucnccd 



Numerically, the HTL approximation gives results somewhat above the 
strictly perturbative result to order g 2 , and, using the same gap equations 
for the asymptotic thermal masses of quarks as in the previous application, 
the result is halfway towards the ideal-gas result, whereas a strictly per- 
turbative inclusion of the plasmon effect would have resulted in a result in 
excess of the free one. 

Fig.ll9lshows the results for quenched QCD together with recent results 
for two different conti nuum extrapolation J * (both are higher than the 
previous lattice resultff^l) , Fig. [30] shows o ur pr edictiorff^H for Nf = 2 
in comparison with the lattice data of Ref. ' \ for which a continuum 
extrapolation is not yet available, however. 

It would be interesting to compare our predictions with HTLPT to two- 
loop order; the one-loop HTL-resummed results from charge correla tor J^ ^ 
suffer from severe overcounting while missing out the plasmon terrr^EII so 
that one should not perform a comparison yet. 



4.5.2. Off-diagonal quark number susceptibilities 

Since lattice results are starting to become available also for the off-diagonal 
QNS, which in the limit of massless QCD we denote collectively by x, we 
briefly discuss this quantity, too. In contrast to the diagonal QNS, the ideal- 
gas value of off-diagonal QNS is zero, so this directly probes the interactions 
in the system. 

The leading-order contribution to off-diagonal QNS require two fermion 
loops connected by gluon lines. The diagram with just one gluon exchange 
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4.5 5 

T/T c 

Fig. 20. Comparison of our results for x/xo m massless Nf=2 QCD with the lattice 
results of Ref.Effl (no continuum extrapolation). Full lines refer to the HTL approxima- 
tion, dash-dotted lines to NLA, when the MS scale is varied from p, = irT to p = AirT. 



vanishes by colour neutrality. The one with two gluon exchange is non-zero, 
but because the fermion loops are then even functions of fi, it contributes 
to x on ly when // ^ 0, starting at order g 3 , namely 



Xij 



g 4 (jV 2 _ 1)Tflifij 

167r 5 m£> 



for i =/= j . 



(112) 



This contribution is indeed contained in the next-to-leading approximation 
discussed above, but only after the corrections to the asymptotic thermal 
masses have been included. 

When all chemical potentials vanish (which is the case most easily ac- 
cessible to lattice studies), the above approximations of the QNS become 
strictly diagonal. In this limit the lowest-order diagram contributing topff- 
diagonal QNS is the 4-loop "bugblatter" diagram shown in Fig. 1211 -^^ 1 In 
this diagram e ach q uark loop produces an effective C-odd vertex for elec- 
trostatic gluon J I which in tu rn gi ve rise to a logarithmically enhanced 



:acn q 

contribution at order g 6 , readin; 



X_ 

Xo 



(N 2 - 1){N 2 -4) (g\G 



(^) 6 logi. 
\7T/ q 



128iV 2 \tt/ - g 

This vanishes in SU(2) gauge theory (though not in , and in ( 

less) QCD this leads to 

5 /o\ 6 , 1 10 



X_ 

Xo 



N=3 



144 



m 6 io g i 

\ir / q 



9tt 3 



a,log(l/a fl ). 



(113) 



mass- 



(114) 
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Fig. 21. Lowest-order diagram in the thermodynamic potential that contributes to off- 
diagonal susceptibilities cWi /d/ij at fi = 0. 

Assuming that the undetermined (though perturbatively computable) 
constant under the log is of order one, the natural order of magnitude of 
y/Y n = y/T 2 is ~ 1CT 4 for T ~ 3T C . Recent lattice studies in N f = 2 
QQL^nni321h a ve f° un( i or ^y va hies consistent with zero so far, with errors 
claimed to be below 10 -6 , but this perhaps remains to be confirmed by 
independent groups before it can be viewed as a new puzzle arising from 
the comparison of lattice and perturbation theory. 

5. Conclusion and perspectives 

In this paper, we have reviewed progress made recently in the calculation of 
the high temperature quark-gluon plasma using weak coupling techniques. 
Such calculations meet with difficulties, many of which are not typical of 
gauge theories, but also occur in scalar field theories. In fact, what deter- 
mine the validity of a weak coupling expansion in thermal field theories is 
not only the strength of the coupling, but also the magnitude of the thermal 
fluctuations which varies according to the relevant momentum scales. 

The thermodynamical functions are dominated by hard degrees of free- 
dom, carrying momenta of the order of the temperature T. Soft modes, with 
typical momenta of order gT, contribute a correction, which however turns 
out to be difficult to calculate, for essentially two reasons. First, the dy- 
namics of the soft modes receive contributions from the hard modes which 
can only be taken into account through all order resummations. Second, 
although the interactions among the soft modes can be treated perturba- 
tively, the corresponding perturbation theory is an expansion in powers of 
g rather than g 2 , as it is for the hard modes. The latter feature is what 
causes the apparent poor convergence of perturbation theory which has 
been observed in calculations up to order g 5 , both in the scalar case and in 
QCD. 

In QCD, because of the unscreened magnetic fluctuations, perturbation 
theory eventually breaks down. But this occurs at order g 6 only, and ac- 
cordingly one could expect the contribution of the magnetic fluctuations to 
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the pressure to be numerically small. This is in fact suggested by detailed 
calculations based on dimensional reduction . Therefore most of this re- 
view has concentrated on the difficulties which are common to all thermal 
field theories, and which involve the momentum scale k ~ gT. 

The various resummation techniques which have been proposed to over- 
come the difficulties of perturbation theory rely to some degree on a separa- 
tion of scale between hard and soft excitations, and follow the general pat- 
tern of effective theories. This is obviously the case for the approaches based 
on dimensional reduction, but is also true of screened perturbation theory 
and ^-derivable approximations, although somewhat more indirectly. In 
dealing with the soft sector, most resummations use the hard thermal loops 
as common building blocks. Non-linear "gap equations" for the calculation 
of "thermal mass" parameters emerge in various forms, and the difficulties 
associated with the ultraviolet renormalizations are dealt with in different 
fashions, depending on the specific approximations. 

In "Screened Perturbation Theory" , and its extension to QCD referred 
to as "HTL-perturbation theory" , the use of the HTL Lagrangian as a basis 
for the resummation has the advantage to preserve manifest gauge symme- 
try, but also the drawback to lead to rather complex calculations beyond 
the lowest order. Besides the HTLs are used for all momenta, soft and 
hard, while they are accurate only at soft momenta. At any finite order in 
this scheme, one meets with temperature-dependent UV divergences whose 
rcnormalization requires the use of "thermal counterterms" . Although they 
get in principle corrected gradually as one includes higher and higher or- 
ders, these features remain unsatisfactory in the low order calculations that 
one is able to perform. 

At least at the level of the formalism, ^-derivable approximations pro- 
vides a conceptually clearer perspective on the general problem of propaga- 
tor renormalization than the simple "add and subtract" approach of SPT. 
The overcounting is automatically avoided by using skeleton expansions, 
and the gap equations emerge naturally, via the variational principle for 
the thermodynamic potential. The HTLs appear as lowest-order approx- 
imations to the solutions of the gap equations, and are used only in the 
kinematical domains where they are legitimate approximations. However, 
the practical implementation of $ derivable approximations is not exempt 
of difficulties: In general, the gap equations are non-local, which compli- 
cates both their UV renormalization, and their numerical resolution. A fur- 
ther complication comes from the fact that the ^-derivable approximations 
violate gauge symmetry. (However, this violation is parametrically under 
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control, which leaves the possibility to construct approximate solutions to 
the gap equations which are gauge-invariant.) 

One of the most obvious advantages of ^-derivable approximations re- 
mains the simplifications they allow in the calculation of the entropy: The 
two-loop calculation for the pressure turns into an effectively one loop cal- 
culation for the entropy, which furthermore is finite once the self-energy 
involved in the calculation has been made finite. The results reported in 
have been obtained however after additional approximations. 
In the scalar 2-loop model, which is exactly solvable, it has been possi- 
ble to test these approximations by comparing against the respective exact 
results. In QCD, where such a comparison is not possible, further refine- 
ments should be studied, in particular to deal with non-local effects in 
the gap equations. Also, the calculation of the entropy is built on real-time 
information, namely the dispersion relation for the hard quasiparticles. Un- 
derstanding the relation with purely Euclidean formalisms, such as that of 
dimensional reduction where this information does not enter in an explicit 
fashion, is certainly worth further investigations. 

In closing, one should perhaps emphasize that most of the calculations 
reported in this review point to the validity of the quasiparticle picture 
of the quark-gluon plasma at high temperature. That is, they support the 
idea that, by suitably dressing the elementary degrees of freedom, one could 
minimize the effects of their residual interactions. This is particularly trans- 
parent in the entropy calculation. The results obtained in this manner for 
the thermodynamic potentials and also for quark susceptibilities seem to en- 
courage this physical picture, which gives hope for the eventual applicability 
of correspondingly improved perturbation theory to dynamical quantities. 
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